Astronomy & Astrophysics manuscript no. final 


©ESO 2013 


March 4,2013 





Nitrogen line spectroscopy in O-stars 
III. The earliest O-stars * ** 

J.G. Rivero Gonzalez 1 , J. Puis 1 , P. Massey 2 ***, and F. Najarro 3 

1 Universitatssternwarte Munchen, Scheinerstr. 1, 81679 Miinchen, Germany, e-mail: jorge@usm.uni-muenchen.de 

2 Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, AZ 86001, USA 

3 Centra de Astrobiologfa, (CSIC-INTA), Ctra. Torrejon a Ajalvir km 4, 28850 Torrejon de Ardoz, Spain 
Received; Accepted 

ABSTRACT 

Context. The classification scheme proposed by Walborn et al. (2002, AJ, 123, 2754), based primarily on the relative strengths of the 
Niv/14058 and Nni/14640 emission lines, has been used in a variety of studies to spectroscopically classify early O-type stars. Owing 
to the lack of a solid theoretical basis, this scheme has not yet been universally accepted though. 

Aims. We provide first theoretical predictions for the N iv/14058/N ni/14640 emission line ratio in dependence of various parameters, 
, and confront these predictions with results from the analysis of a sample of early-type LMC/SMC O-stars. 

. Methods. Stellar and wind parameters of our sample stars are determined by line profile fitting of hydrogen, helium and nitrogen 

lines, exploiting the helium and nitrogen ionization balance. Corresponding synthetic spectra are calculated by means of the NLTE 
atmosphere/spectrum synthesis code fastwind. 

Results. Though there is a monotonic relationship between the N iv/N m emission line ratio and the effective temperature, all other 
parameters being equal, theoretical predictions indicate additional dependencies on surface gravity, mass-loss, metallicity, and, par- 
ticularly, nitrogen abundance. For a given line ratio (i.e., spectral type), more enriched objects should be typically hotter. These basic 
' predictions are confirmed by results from the alternative model atmosphere code cmfgen. 

Cln, The effective temperatures for the earliest O-stars, inferred from the nitrogen ionization balance, are partly considerably hotter than 

indicated by previous studies. Consistent with earlier results, effective temperatures increase from supergiants to dwarfs for all spectral 
types in the LMC. The relation between observed N iV/14058/N ni/14640 emission line ratio and effective temperature, for a given 
luminosity class, turned out to be quite monotonic for our sample stars, and to be fairly consistent with our model predictions. The 
scatter within a spectral sub-type is mainly produced by abundance effects. 

Conclusions. Our findings suggest that the Walborn et al. (2002) classification scheme is able to provide a meaningful relation 
between spectral type and effective temperature, as long as it is possible to discriminate for the luminosity class. In terms of spectral 
morphology, this might be difficult to achieve in low-Z environments such as the SMC, owing to rather low wind-strengths. According 
to our predictions, the major bias of the classification scheme is due to nitrogen content, and the overall spectral type-T eff relation for 
■^j- , low-metallicity (e.g., SMC) O-stars might be non-monotonic around 03.5/04. 
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J^D ' 1. Introduction their surroundings, e.g., their input of (ionizing) radiation, wind- 

Z^l \ energy and momentum, and nuclear processed material. 

Though important for the evolution of the early Universe and Particularly uncertain is the situation for the earliest O-stars 

the present cosmos, massive stars are not as thoroughly under- (earlier than Q ^ see below)> both with respect to their physical 

stood as desirable to safely infer or predict their interaction with parame ters, and the relation between their spectroscopic defini- 

' tion and these parameters. Since the earliest O-type stars include 

s_i j. „ ~, ', r . , _ , the most massive stars in the Universe, i.e., the top end of the 

3 Based on (l) observations collected at the European Southern „ . . . , n . , , . r . 

Observatory Very Large Telescope, under programmes 68.D-0369, Stellar lnltlal maSS function, such a lack ofknowledge is intol- 

171.D-0237 (FLAMES), and 67.D-0238, 70.D-0164, 074.D-0109 erable. As pointed out by |Massey et al.| d2005|), already a 10% 

(UVES); (ii) observations made with the NASA/ESA Hubble Space uncertainty in the effective temperature, 7/ eff (which is still pos- 

Telescope, obtained from the Data Archive at the Space Telescope sible for the hottest stars to present date), results in a factor of 

Science Institute, which is operated by the Association of Universities two or more uncertainty in the Lyman flux, affecting our under- 

for Research in Astronomy, Inc., under NASA contract NAS 5-26555. standing of the ionization balance of H n regions and, e.g., the 

These observations are associated with programmes 6417, 7739, and porosity o f the interstellar medium (e.g JOev & Kennicutfl 19971 

9412; and (iii) observations gathered with the 6.5m Magellan telescopes Q ey 2006) 



at the Las Campanas Observatory, Chile. 



This large uncertainty is produced because the standard ap- 



** Appendix A and B are only available in electronic form at , , , ,,,,,, ... 

irn — a — ■ 1 proach to derive T e fr, exploiting the Hei/Hen ionization equi- 

http://www.edpsciences.org £ eIt ' K „ 6 " _ . . , M . 

*** Visiting astronomer, Cerro Tololo International Observatory, llbrlum ' falls ln . the earllest " star re 8 lme ' Thls 18 the result 

National Optical Astronomy Observatory, which is operated by the of rather insensitive He n lines to changes in r eff , and quite 

Association of Universities for Research in Astronomy (AURA) under weak strategic He I lines, with their equivalent widths falling 

cooperative agreement with the National Science Foundation. below 200 mA, about the limit achievable with normal photo- 
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graphic emul sions in the 197 0s. The apparent absence of Hei 
was used by Walborn| (Il97ll) to extend previous classification 
criteria to the 03 spectral type, which then displayed a degener- 
acy with respect to effective temperat ure. This degeneracy was 
partially broken by iKudritzkil d 19801) and ISimon et all d 1983b 
who used better photographic data to detect He i in the spec- 
tra of some 03 stars, but even today with modern CCDs and 
high S/N the detectio n of He i is not alway s achievable. To cir- 
cumvent this problem Walbor n et al.l (|2002) suggested to use the 
N iv /14058/Nm^46403 (hereafter Nrv/Nm) emission line ratio 
as the primary classification criterion for the earliest spectral 
types, instead of the He i/He n absorption line ratio. By means of 
this criterion, the former 03 class was split into three different 
types 02, 03, and 03.5, being 02 the degenerate one instead. 



This scheme has been used in a variety of studies during the 
past years to classify spectra of early O-stars in the Large and 
Smal l Magellanic Clouds ( LMC and SMC) (e.g..lWalborn et alj 
12001 IMassev et all l200l 120051 [20091 lEvans et alJ 120061) and 



in the Milky-Way, e.g.,[Sotaet al. (2011). However, these addi- 
tional sub-types are not ( yet) universally accepted. In particular 
IMassev et al.l d2004l 120051) criticized that relying on the relative 
strengths of the optical nitrogen emission lines lacks a solid the- 
oretical basis, since the (photospheric) line emission is the result 
of complex NLTE processes. 



Table 1. Model-grid used within this work: coverage of funda- 
mental parameters. For T e g and log g coverage, see also asterisks 
in Fig.|S] 



Parameter 


Range 


Typical step size 


r cff (kK) 


35.0- 55.0 


1.0 


logg (cgs) 


3.2,3.5-4.2, 4.5 


0.1 


log Q" /series 


-14.00/A, -13.50/B, -13.15/C 


0.35 




-12.80/D, -12.45/E, -12.10/F 






0.08, 0.10-0.20 


0.05 


Galactic 


Z 


z Q 




[N] 


7.64, 7.78 - 8.98 


0.2 


LMC 


Z 


0.5Z G 




[N] 


6.90, 6.98 - 8.58 


0.2 


SMC 


Z 


0.2Z o 




[N] 


6.50, 6.78 - 8.38 


0.2 



<n) in units of M yr-7(km s- l R Q ) [ 5 . 

Notes. Other model parameters adopted as follows: wind terminal ve- 
locity, Vqq, as a functi on of the photospheric escape velocity, v esc (see 
iKudritzki & Puisl l2000l) : stellar radius, R,, as a function of spectral 
type and luminosity class, corresponding to prototypical values; veloc- 
ity field exponent, fj = 0.8; and micro- turbulence, v m j c = 10 kms -1 . 
Nitrogen b aselines abundances (first entry in [N] range) have been 
drawn from Hunter et al. (20071). 



To provide more insight into this and related matter, we 
started a series of publications dealing with nitrogen spec- 
troscopy in O-type stars. I n the first paper of this series 
dRivero Gonzalez et al.l 1201 ll hereafter Paper I), we concen- 
trated on the formation of the emission at N m/1/14634 - 4 640 - 
4642. In the follow-up paper (Hi vero Gonzalez et al II20T2L here- 
after Paper II), we investigated the N iv/14058 emission line for- 
mation, and applied our accumulated knowledge to derive nitro- 
gen abundances for an LMC O-star sample. 

The primary goal of the present paper is a quantitative study 
of the atmospheric parameters of the earliest O-stars, by means 
of nitrogen line spectroscopy and building on the results from 
our previous work within this ser i es. P articularly, we concen- 
trate on testing the Wal born et alj (|2002|) classification scheme 
on its capability of providing a reasonable relation between spec- 
tral types and effective temperatures. To this end, we investigate 
the theoretical N iv/N m emission line ratio, and the impact of 
different parameters on this ratio. Subsequently, our predictions 
are confronted with corresponding observational results, derived 
from an analysis of an early O-type sample of LMC/SMC stars. 

This paper is organized as follows. In Sect. [2] we describe the 
tools used within this work, the atmospheric code fastwind and 
a suitable model grid. Open questions from our previous studies 
regarding the formation of the N m and N iv emission lines are 
addressed in Sect. [3] Section [4] presents first theoretical predic- 
tions on the N iv/N m emission line ratio. In Sect. [5] we compare 
our results and predictions with cor responding ones from t he al- 
ternative atmospheric code cmfgen dHillier & Miller 1998). The 
stellar sample and the observations used within this study as well 
as the procedure to determine stellar/wind parameters together 
with nitrogen abundances are presented in Sect. |6] Finally, we 
provide a discussion of our results in Sect. [7] and summarize our 
findings and conclusions in Sect. [8] 



2. Code and model grid 

All calculations within this work were performed by means 
of the NLTE atmosphere / spectrum synthe sis code fastwind 
(ISantolava-Rev et aLII 19971 iPuls et akl l2005). using the recently 
updated version vlO. 1 (see Paper II). Most results presented here 
(except for the fine-tuned fits in Sect. [6} are based on a model- 
grid, with H, He, and N as 'explicit' elements. Corr esponding 
model atoms have been described in IPuls et al" I (12001 (H/He) 
and Papers I/II (N), and Table Q] provides the coverage of im- 
portant grid parameters. Details of the basic grid were already 
provided in Paper II, and only its 'hot' range (T e g > 35kK) was 
used for our current analysis. This subgrid has been extended 
to cover a broader range in background metallicity, Z = 1, 0.5, 
0.2 Z , associated with the Milky Way (MW), Large Magellanic 
Cloud (LMC), and Small Magellanic Cloud (SMC), respec- 
tively. Moreover, we increased the sampling with respect to 
log g, and also the coverage of the wind-strength parameter (or 
optical depth invariant), Q = M/iy^R*) 1 - 5 , towards larger values, 
resulting in 6 model series of different wind-strength (from se- 
ries 'A' with log Q = -14.00 to series 'F' with log Q = -12.10, for 
details see Table[T]l. Finally, the sampling in helium content, Yu e 
= Nhs/Nh, and nitrogen abundance, [N] = log l0 (N^/NH) + 12, 
has been improved, for studying the reaction of important nitro- 
gen lines on extreme changes in [N] (Sects. 1331 and l4~2~3l l. and 
for better constraining abundances in our analysis of LMC/SMC 
stars (Sect. [6j. The current grid accounts for a total of 104,000 
models. 

The major potential of such a model-grid (besides theoretical 
investigations) is the possibility of obtaining rather precise esti- 
mates of stellar parameters within a reasonable amount of time 
(few minutes). With the advent of large stell ar surveys, such as 
the VLT-FLAMES Tara ntula survey (VTFS, lEyans et al.ll201 ll) 
or the IACOB project (Sim on-Diaz et al.l [20 1 ll) . this will turn 

1 The individual components NIILI/14640.64-4641.85 become 
blended for a projected rotational velocity v sin i > 40 km s . 
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H_eps 



H_delta 




3960.0 3966.7 3973.3 3980.0 
wavelength (A) 



H -alpha 



E 0.8 




6543 6S58 6570 6583 
wavelength (A) 

HeI4713 




4703.0 4709.7 4716.3 4723.0 
wavelength (A) 

HeII6406 




E 0.96 



6396.0 6402.7 6409.3 6416.0 
wavelength (A) 

NIII4097 




4096.3 4096.6 4098.0 4099.3 
wavelength (A) 

N1I14510-14-18 





H _gamma 



H_beta 



4090.0 4096.7 4103.3 4110.0 
wavelength (A) 

He!+He!I4026 




4016.0 4022.7 4029.3 4036.0 
wavelength (A) 

HeII4200 




).0 4196.7 4203.3 4210.0 
wavelength (A) 

HelI6527 




6517.0 6523.7 6530.3 6537.0 
wavelength (A) 

NIII4195 




4193 4195 4197 4199 
wavelength (A) 

N1V4058 





4330.0 4336.7 4343.3 4350.0 
wavelength (A) 

HeI4387 




4377.0 4383.7 4390.3 4397.0 
wavelength (A) 

HeII4541 




4531.0 4537.7 4544.3 4551.0 
wavelength (A) 

HeI6678+HeII6683 




6673.0 6679.7 6686.3 6693.0 
wavelength (A) 



NIII4379 



1.02 



1.01 



1.00 



g 0.99 



0.98 




4376 4378 4380 4382 
wavelength (A) 

NIV6380 




4508 4512 4518 4520 
wavelength (A) 



4051.0 4055.7 4060.3 4065.0 
wavelength (A) 



0.90 



6373.0 8378.3 6383.7 6389.0 
wavelength (A) 




4860.0 4858.7 4863.3 4870.0 
wavelength (A) 

HeI4471 




4461.0 4467.7 4474.3 4481.0 
wavelength (A) 

HeII4686 




4676.0 4682.7 4689.3 4696.0 
wavelength (A) 

NIII4003 




4000 4002 4004 4006 
wavelength (A) 

NIII4634-4640-4642 





1.08 




1.06 


i 




Q. 


1.04 






C 




V 

2> 


1.02 


V 




S 

V 


1.00 




0.98 



4828 4834 4640 4646 
wavelength (A) 

NV4603-4619 




S 0.96 



0.94 



4597 4806 4616 4625 
wavelength (A) 



Fig. 1. High resolution optical spectrum of the Galactic 03 V((f*)) star HP 64568 (green), compared with the best-fitting synthetic 
H/He/N spectrum from our grid (black, convolved with v sin i = 100 km s , see lMarkova et al EoTTl) . N iv^6380 and N v.14603 not 
observed. Grid parameters: T e s = 48 kK, logg = 4.0, log Q = -12.8 (series 'D'), iHe=0.1, [N] = 8.38. Fine tuning of the parameters 
can improve the fit. 
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out to become a very useful tool. As an example for the present 
quality of our models, we show in Fig. [T]the comparison of a 
high resolution (resolving power R - 48,000), high S/N (> 150) 
spe ctrum of the Galactic 03 V((f*)) star HD 64568 (for details, 
see lMarkova et al.ll201 lb with the best fitting synthetic spectrum 
from our grid models. Obviously, the grid resolution is sufficient 
to achieve quite a fair representation of the observed spectrum, 
both with respect to H/He and N (note that lines from N m, N iv, 
and N v are present in parallel). Further analyses of Galactic ob- 
jects will be performed in a forthcoming paper. 

3. Nm/Niv emission lines - parameter-dependence 

Our previous studies on the formation of N mvl/14634 - 4640 - 
4642 and N iv/i4058 prompted a number of questions, which are 
investigated in the following. First we concentrate on the Nm 
triplet emission in low-Z environments, to discriminate the coun- 
teracting effects of line blocking (less blocking - more emis- 
sion) and mass-loss (less mass-loss - less emission) under these 
conditions. Subsequently, we study the impact of Z and [N] on 
NivA4O580 

3.1. Nm emission line formation: EUV line-blocking vs. 
wind-strength 

In Paper I we argued that, for Galactic conditions, the 
canonical explanation for the presence of emission at 
Nni/ 114634 - 4640 - 4642 (r elated to dielectronic recombina- 
tion, |MSalas^&Hummea[l973|) no longer or only partly applies 
if one accounts for the presence of line-blocking/blanketing and 
winds. The key role is now played by the stellar wind, as long 
as the wind-strength is large enough to enable a significantly ac- 
celerating velocity field already in the photospheric formation 
region of the resonance line(s) connected to the upper level of 
the involved transition. Furthermore, our study implied that par- 
ticularly the efficiency of the 'two electron' drain (depopulating 
the lower level, see Paper I) is strongly dependent on the de- 
gree of EUV line-blocking, i.e., on Z. For a given wind-strength 
and nitrogen abundance, the emission should become stronger 
in low-Z environments, because of less blocking (see Fig 16 
in Paper I). Nevertheless, this comparison might be unrealis- 
tic, since less blocking goes hand in hand with a lower wind- 
strength, which might (over-) compensate the discussed effect. 
Moreover, one might need to consider a lower nitrogen content, 
owing to a lower baseline abundance. 

Here we investigate the combined effect. First, we compare 
the behavior of the N in emission lines for two different Z en- 
vironments (MW and SMC), applying a consis tent scaling of 
log Q, via M oc (Z/Z Q ) - 72 (clumping corrected, iMokiem et al.l 
2007b) and Vco oc (Z/Z ) 013 (iLeitherer et al.ll992l) . For the SMC 
with Z/Z = 0.2, this yields a reduction of log Q by -0.37 dex 
compared to the Galactic case, and corresponds well to the 
step size of A log Q = 0.35 used within our model grid. Thus, 
Galactic models from, e.g., series 'E' need to be compared with 
SMC models from series 'D'. 

In Fig. |2] (upper panel), we compare the equivalent widths 
of Nmvl/14640.64, 4641.853, in the following abbreviated by 
N ni/14640, as a function of T e ff, for MW series 'E' (black/solid) 
and SMC series 'D' models (red/dashed). All models have 
the same gravity, logg = 4.0, and the same nitrogen content, 

2 For a detailed description of the formation process of N m and N iv 
emission lines, we refer to Paper I and Paper II, respectively. 

3 The blue component of the triplet shows a similar behavior. 
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-0.3 L , . . . . , , . . , , . , . . , , . . . : 

3.5x10" 4.0x10" 4.5x10" 5.0x10" 5.5x10" 

T.„ (K) 

Fig. 2. Equivalent width of Nni/14640 as a function of T e s, 
for MW and SMC models at logg= 4.0 and mass-loss rates 
scaled according to Z (see text). Solid/black and dashed(- 
dotted)/red curves refer to MW series 'E' and SMC series 
'D' models, respectively (A log Q = -0.35 dex). Upper panel: 
[N] = 7.78 for both Galactic and SMC models. Middle panel: 
comparison accounting for the theoretically expected maximum 
r N 1 enrichment, draw n from tailored evolutionary calculations 
bv lBrott etal] d2011al) : MW models with [N] =8.18 and SMC 
models with [N] = 7.58. Lower panel: SMC series 'D' models, 
with [N] = 7.58 (dashed-dotted) and [N] = 7.78 (dashed). The 
line emission increases with increasing [N]. 

[N] = 7.78 dex (almost solar, lAsplund etal1l2005l l2009h . As 
usual, we define equivalent widths (hereafter EWs) to be pos- 
itive for absorption and to be negative for emission lines. The 
result of this comparison is similar to our findings from Paper I. 
Even when we account for consistently scaled mass-loss rates, 
the low-Z models result in more emission. Thus, lower mass- 
loss rates associated with low-Z environments do not compen- 
sate the increase in emission strength due to a lower degree of 
line-blocking. 

So far, we neglected the fact that a lower Z also implies 
a difference in the [N] baseline abundance, which needs to be 
considered in a final comparison. In the middle panel of Fig. [2] 
we compare the same series of models as in the upper one, ac- 
counting now for more consistent abundances (MW: black/solid, 
SMC: red/dashed). Here we used [N] values in agreement with 
the theoretically ex pected maximum enrichment as provided by 
iBrottet al.l(l2011al) . forMW and SMC models of 40 M and ini- 
tial rotation velocities of 270 kms -1 . This roughly corresponds 
to models with [N] = 8.2 and 7.6 for the MW and SMC, respec- 
tively. Q 

4 Unfortunately, there are only few abundances studies in the O-star 
regime to confirm these assumptions , and most of them are biased to- 
wards late O-types (reviewed by, e.g. Herrero 2003, Herrer o~& Lennor] 
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3.5x10" 4.0x10" 4.5x10" 5.0x10" 5.5x10" 

T.„ (K) 

Fig. 3. Equivalent width of Nni/14640 (upper panel) and 
NiV/14058 (lower panel) as a function of T eS , for MW (black) 
and SMC (red) models, at logg=4.0 and [N]=7.78 dex. Solid 
lines: series 'A' (thin winds), dashed lines: series 'E' (prototypi- 
cal for Galactic supergiants). 



It turns out that with the inclusion of more realistic abun- 
dance conditions, the increase of emission strength due to less 
blocking becomes strongly suppressed, and now the MW mod- 
els display stronger emission that the SMC ones. This might ex- 
plain the relatively low number of "f" objects among the SMC 
O-stars (see Sect. 17.6b . 

For convenience, the lower panel of Fig. [2] displays a di- 
rect comparison of SMC models with expected maximum en- 
richment ([N] = 7.58, dashed-dotted) and a solar nitrogen con- 
tent ([N] = 7.78, dashed). Obviously, the emission strength of 
N ni/14640 increases with increasing nitrogen content (see also 
Fig. 15 in Paper I). 

3.2. A/iv/14058 - dependence on background abundance 

Though we studied the response of the N in triplet on different 
background abundances already in Paper I, a similar analysis 
for NiV/14058 is still missing. This is now done in Fig. [3] by 
means of our model-grid. We display models with thin winds 
(series 'A' ; solid) and with wind-strengths typical for Galactic 
supergiants (series 'E'; dashed), and compare the impact of MW 
(black) and SMC (red) background abundances. All models have 
the same gravity, logg= 4.0, and the same [N] = 7.78. To en- 
able a better comparison, we show the effects for both N ni/i4640 
(upper panel) and N iv/14058 (lower panel). We clearly see that 
N ni/14640 is much more influenced by Z than NiV/14058, irre- 
spective of wind-strength and r e jf . For LMC background abun- 
dances, Z = 0.5 (not displayed), the effect on N iv/14058 becomes 
almost negligible. 

With respect to the formation of this line, all A' and also 
the 'cooler' 'E' models (with r e ff < 45 kK) behave as discussed 
in paper II. A lower background Z induces a more depopulated 
ground state, owing to somewhat higher ionizing fluxes around 

2004, iMorell 120 09). Hea p et all j2006ft . following previous work 
by iBouret et al.l d2003l> . found most of their 18 SMC sample stars 
to be enriched, and more than hal f of the sample disp layed [N] 
> 7.5. Regarding G alactic objects, Mar tins et al.l d2012bh (see also 
Martins et al. 2012a) recently found a typical enrichment by 0.4-0.6 dex 
above the Galactic baseline abundance, corresponding to [N] x 8.0-8.2. 
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-0.6 I , i , , , , i , , , , i , , , , i , , , , I 

3.5x10" 4.0x10" 4.5x10" 5.0x10" 5.5x10" 

T.„ (K) 

Fig. 4. Equivalent width of Niv/14058 as a function of T e ff for 
models with [N] = 7.78 (black) and [N] = 9.0 (red). All mod- 
els correspond to log,g=4.0 and Z = Z Q , with mass-loss rates 
according to series A' (solid) and 'E' (dashed). See text. 

the N iv edge at A — 160 A. Consequently, the drain of the lower 
level of the /14058 transition, 3p, via the 'two electron' tran- 
sitions becomes enhanced, and more emission is produced at 
lower Z. For 'E' models with r e ff > 45 kK the situation changes, 
and now the higher Z (MW) models produce slightly more emis- 
sion: At higher r e jf and M, the resonance line towards 3p (at 
A = 247 A) leaves detailed balance and 3p becomes pumped, 
stronger at low Z because of higher EUV fluxes. Thus, less emis- 
sion at NiV/14058 is produced, compared to a Galactic environ- 
ment. 

Overall, however, background abundances have a weak ef- 
fect on N iv/14058, much weaker than wind-strength effects (e.g., 
red solid vs. red dashed). This is related to the much higher sen- 
sitivity of the Niv continuum on mass-loss rate. Accordingly, 
also the helium content has a certain impact on the N iv emis- 
sion strength, since this parameter controls the overall flux level 
of the Hen continuum including the Niv edge. In particular, an 
increase of Yn e decreases the N iv emission strength. 

As already discussed above and in Paper I, N ni/14640 reacts 
much stronger on different background abundances (e.g., upper 
panel of Fig. [3] black vs. red ), and in a similar way for all mass- 
loss rates. Because of this different behavior, the (theoretical) 
NiV/14058/N ni/14640 emission line ratio, studied in Sect. [4] is 
affected by metallicity. 

Interestingly, the reaction of N ni/14640 on M becomes neg- 
ligible for models with T e ff > 46 kK. At these temperatures, N m 
has become a real trace ion, and the (weak) line emission is 
formed in quite deep layers, hardly affected by mass loss and 
velocity field (see also Fig. 10 in Paper I). In this case, the rel- 
ative overpopulation is caused by recombination cascades and 
depopulation of the lower level by 'two electron' drain (partic- 
ularly for low Z conditions), whilst dielectronic recombination 
remains unimportant. 

3.3. Niv/14058 - dependence on nitrogen abundance 

Figure|4]displays the reaction of N iv/14058 on variations of [N] 
and log Q for MW models. As in Fig. [3] solid and dashed lines 

5 A similar behavior is found for the N in triplet emission, for T eff < 
40 kK. 
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Fig. 5. Radiation temperatures and NLTE departure coefficients for two MW-models at r e jf= 45 kK, log,g= 4.0, and a low density 
wind (series 'A'), with solar (black) and strongly enhanced nitrogen ([N] = 9.0, red). Left panel: Enhanced nitrogen leads to lower 
Niv (and Nm) continuum-fluxes, here expressed as radiation temperatures. Corresponding edges at 160 A and 261 A indicated by 
vertical lines. Right panel: The lower fluxes give rise to less depopulated N iv ground- (solid) and 2p 2 (dotted) states, where the 
latter are the lower levels of the important two-electron transitions draining N iv 3p. Consequently, there is more absorption/less 
emission at Niv/14058 when [N] becomes increased. Formation region of N iv/14058 indicated by dashed lines. 



correspond to model series A' and 'E', respectively. We com- 
pare models with [N] = 7.78 (black) and highly enhanced nitro- 
gen, [N] = 9.0 (red), selected to demonstrate extreme effects. 

First, we concentrate on the influence of [N] for model series 
'E' (supergiant mass-loss rates, dashed). These models behave as 
expected. As for the N m triplet (e.g., Fig. [2] lower panel), we ob- 
tain more emission when we increase the nitrogen content. This 
is mostly because the formation zone, due to the increased num- 
ber of absorbers/emitters, is shifted outwards into the transition 
region photosphere/wind where the relative overpopulation be- 
comes very large or even inverted (cf. Fig. 4 in Paper II). Small 
inaccuracies in the population ratio (e.g., due to inappropriate 
gridding) can lead to sizeable effects in line-strength when close 
to inversion, and this is the reason for the non-monotonicity in 
EW encountered for the high [N] 'E'-models. 

Somewhat unexpectedly, however, we found the opposite 
reaction for low-M models (solid). For large nitrogen content 
(red), we either obtain more absorption or less emission than 
solar- [N] models (black), for almost the whole temperature 
range. Furthermore, the turning point from absorption to emis- 
sion occurs at hotter T e g (by 5 kK) than for the solar- [N] models, 
whilst for the 'E' series this turning point remains rather unaf- 
fected by nitrogen content. The underlying mechanism is again 
related to a diminished drain of the lower level when [N] be- 
comes increased. Here, a higher [N] is responsible for a lower 
Niv continuum-flux (Fig. [5] left), leading to less depopulated 
ground- and 2p 2 states (Fig. [5] right), and thus to a higher popu- 
lation of 3p, i.e., to more absorption/less emission at N iv/14058. 
Note that this process is similar to the effect produced by a higher 
Z (Sect 153b. 

These findings imply an important consequence. According 
to our predictions, in a certain T e ff range Niv/14058 might 
switch, for growing [N], from emission to absorption, provided 
the wind-strength is not too high! The other way round (and ex- 
ploiting the results from the complete grid): If N iv/14058 is ob- 
served in absorption at 44 kK < T e g < 50 kK (for Galactic stars), 
this would be an indication of strong nitrogen enrichment (and 
mass-loss rates below 'D' corresponding to log Q < -12.8). The 
lower limit in T e g corresponds to the absorption/emission turn- 



ing point for 'D' models (at higher T e ff and similar or higher M, 
only emission lines are predicted, independent of [N]), whilst the 
upper one refers to the same point for model series A' . 

4. Predictions on the Niv/Nm emission line ratio 

4.1. Overview 

The complete Wal born et alJ (I2002I) classification scheme pro- 
posed for 'normal' O-starsf] which also covers the 04 type, 
is summarized in Table [2] Walborn et al. suggested to use the 
Niv/Nm emission line ratio as the primary classification cri- 
terion for the earliest spectral types, instead of the He i/He n 
absorption line ratio. This scheme has been used in a variety 
of studies during the past years to classify spectra of early O- 
stars (e.g.JWalborn et al.ll2004l iMassev et al.ll2004t I20051 12009L 
lEvans et aUl2006l ISota et aLll201 lb . though there are still some 
controversial issues, (i) The classification criteria are not quanti- 
tative and involve secondary stat ements regar d ing th e strength 
of He 1/14471, see Table |2] (ii) iMassev et all d2005l hereafter 
Mas05) have criticized that relying on the strength of the nitro- 
gen emission lines lacks a theoretical basis, i.e., it is not clear 
whether T e s is the only parameter that differentiates the newly 
defined spectral types. Indeed, our previous and present work 
implies that the emission strength of N iv/14058 (and also that 
of N ni/14640, at least until T e g » 46 kK) crucially depends on 
M. Based on their analysis of early LMC and SMC stars, Mas05 
pointed out that for stars with similar T e ff and log g the N iv/N m 
ratio could vary by the full range defined for the scheme (but see 
Sect. 16.4b . Furthermore, they suggested that any spectroscopic 
classification should be able to constrain T e g without knowledge 
of other important p arameters such as log g or M. (iii) Already 
Walb orn et al.1 (120021) pointed out that even though there were no 
indications of any correlation between the newly defined spec- 
tral types and the corresponding host galaxy, effects of Z on the 
emission line ratio might need to be considered, accounting for 

6 ICrowther & W alborn ( 201 1 ) have updated the classification scheme 
for 02-3.5 If*/WN5-7 stars using the morphology of H^. 
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Table 2. Classification scheme for spectral types 02-04 using 
the emission line ratio N iv/14058/ N m/14640 a nd the strength of 
He 1,1447 1 . as defined bv lWalborn et all d2002l). 



Spectral type 


Criteria 


Supergiants 


02 If* 


N iv » N m, no or very weak He i 


03 If* 


N iv > N m, very weak or no He i 


03.5 If* 


N iv ~ N m, very weak He i 


04If+ 


N iv < N m, weak He i 


Giants 


02 III(f*) 


N iv » N m, no or very weak He i 


03 m(f*) 


Niv > Nm, very weak or no He i 


03.5 m(f*) 


N iv ~ N m, very weak He i 


04 ni(f+) 


N iv < N m, weak He i 


Dwarfs 


02 V((f*)) 


N iv » N m, no He i 


03 V((f*)) 


N iv > or ~ N m, very weak He i 


03.5 V((f+)) 


N iv < N m, very weak He i 


04 V((f+)) 


no N iv, weak He i 



Notes. Luminosity classes defined as follows. Supergiants (I): 
Heii/t4686 in emission; giants (III): Heii/t4686 in weak absorption/P- 
Cygni profile; dwarfs (V): Heii/14686 in strong absorption. Note that 
the f+ designation recently became obsolete since the Si iv emission at 
,M4089 - 4116 has been establ ished as a common feature in normal 
O-type spectra dSota et al.ll2011h . 

the results from ICrowtherl d200 ) for WR-stars. Using synthetic 
WN models, Crowther had found that earlie r spectral types are 
predicted at lower metallicity, following the ISmith et al.l (fl996) 
classification scheme for WN stars. Since both nitrogen emis- 
sion lines seem to react differently on variations of Z (Sect. l3T2l >. 
it is clear that such a dependence cannot be ruled out. 

In the following, we will address these and other problems, 
first by means of theoretical predictions for the N iv/N m emis- 
sion line ratio, and later on by a comparison with observed spec- 
tra. 

4.2. Basic considerations 

We took advantage of the large model-grid described in Sect. [2] 
and analyzed the influence of various parameters (Z, log Q, [N]) 
by studying the iso-contours of specific emission line ratios in 
the r e ff-logg plane (Figs.[6]to[8]l. Here and in the following, we 
always used models with l^e =0.1. 

To discriminate the specific spectral types by nitrogen lines 
alone, one has basically to account for five different, qualitatively 
defined ranges with respect to the line strengths of N iv/i4058 
vs. N ni/14640 (see Table[2]i. Note that the ranges for luminosity 
classes I/III are somewhat shifted relative to class V. To allow for 
a quantitative description, we investigated the behavior of three 
extreme values, namely Niv/Niii = 0.1 (lower limit for 04 I/III 
and 03.5 V), Niv/Niii = 1 (03.5 I/III and 03 V) and Niv/Niii 
= 10 (representative for 02 I/III/V). 

Before going into the details of our analysis, we compare in 
Fig. [6] the Niv/Nm emission line ratios expressed in terms of 
EW (left) and line-strength (right, quantified in terms of emis- 
sion peak height), to ensure that different definition^ would 
not lead to different conclusions. Except for small subtleties, 
there is no significant change in the run of the iso-contours, 

7 Note that line ratios for He i/He n are partly defined using EWs. 



regardless whether EWs or line-strengths are used to derive 
the line ratio. The encountered subtle differences are mostly 
caused by the considerable wind-strength (model series 'E') 
used in Fig. [6] Particularly for NiV/14058, enhanced EWs ow- 
ing to extended wings are produced, whilst the peak heights re- 
main unaffected. We have also convinced ourselves that a con- 
volution of the theoretical spectra with typical rotational speeds, 
v sin i = 100 kms" 1 , and/or a degrading to a resolving powe r 
of 6000 have a rather low impact (but see iMarkova et al ■1 I20TT1) : 
Generally, the iso-contours are shifted to somewhat lower T e g, 
by roughly 500 to 1000 K in extreme cases of high wind-strength 
and high [N] . In the following we conce ntrate on the ratio of lin e- 
strengths alone, as originally defined by Walborn et aD (l2002l) . 

The inspection of the different iso-contours displayed in the 
r e ff-logg plane, Figs. [6] to [8] allows us to infer an important 
characteristics for the emission line ratio. As it is true for the 
Hei/Hen line ratio (e.g., Mas05), also the Niv/Niii line ratio is 
a sensitive function of surface gravity. The hotter the tempera- 
ture, the larger the log g-value necessary to preserve a specific 
ratio (ionization vs. recombination). This trend seems to have 
vanished for the cooler models (T e g < 40 kK) at high logg in 
Fig. [6] but here the influence of gravity is counteracted by the 
rather strong wind (see Sect. 14.2.2V 

4.2.1. Impact of background metallicity 

Figure [6] shows the dependence of the line ratio on T e ff and logg 
for a fixed wind-strength, log Q = -12.45, and [N] = 7.78, for 
iso-contours corresponding to N iv/N m = 0. 1, 1 and 10. To study 
the influence of Z, we display the predicted behavior for MW 
(red), LMC (green), and SMC (blue) O-stars. 

Even a first inspection indicates the potential of this line 
ratio as a temperature diagnostics. For a fixed surface gravity, 
e.g., logg = 4.0, we find difference on the order of 3 kK (for 
MW models) to 8 kK (for SMC models) between iso-contours at 
N iv/N m = 0. 1 and 1 , roughly corresponding to 04 (lower limit) 
vs. 03. 5 I/III or 03 V. For low-M models this spread becomes 
somewhat smaller, 6 kK for SMC stars, whilst for Galactic ob- 
jects the Niv/Nm = 0.1 iso-contour is no longer present since 
Niv turns into weak absorption (Fig. [3}. The difference between 
the Niv/Nm = 1 and 10 iso-contours (03.5 I/III or 03 V vs. 
02) is on the order of 4-5 kK for all metallicities, decreasing to 
2.5-3 kK for low M. 

The specific differences of the emission line ratios as a 
function of Z can be easily explained in terms of our conclu- 
sions from Sect. l3.2l and by remembering the basic features dis- 
played in Fig. [3] First, N ni/14640 is strongly affected by Z (more 
emission at lower Z, irrespective of temperature), in contrast to 
N iv/14058. The latter line is increasing in strength over almost 
the complete range in T e ff , whilst the former is increasing until 
a certain maximum located at « 40 kK, and decreasing there- 
after. Taken together, this implies that lower are required to 
produce similar line ratios for low-metallicity objects at cooler 
temperatures, compared to Galactic objects (lower T e g means 
less N in emission, to compensate for the basically larger emis- 
sion due to lower Z). Beyond 40 kK, on the other hand, higher 
r e ff are required instead to reduce the increased N in emission. 
This explains why the SMC/LMC iso-contours at N rv/N m = 0. 1 
are located at cooler than their Galactic counterpart, whilst 
the corresponding Niv/Niii = 1 and 10 iso-contours are on the 
hotter side. Because the Z effect is largest for SMC objects, also 
the T e g differences between the N iv/N in = 0.1 and 1 levels are 
largest in this situation. 
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NIV4058/NIII4640(EW): series E, log Q = -12.45 NIV4058/NIII4640(strength): series E, log Q = 
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Fig. 6. Niv/14058/N ni/14640 for Galactic (red), LMC (green) and SMC (blue) O-stars. Displayed are iso-contours of the emission 
line ratio in the T^Aogg plane (model series 'E', [N] — 7.78), for values of Niv/Nm = 0.1 (dashed), 1 (solid) and 10 (long- 
dashed). The asterisks indicate the position of the grid-models used. Left panel: emission line ratio calculated using EWs; right 
panel: emission line ratio calculated using line-strengths. 



Finally, let us mention that the influence associated with Z 
(around 2-3 kK for a given line ratio) is comparable with the 
typical spread present in early spectral types, and thus in agree- 
ment with the non-detection of any correlation between the new 
spectral types and the host galaxy (see above). 

4.2.2. Impact of wind-strength 

To test the influence of log Q on the emission line ratio, in Fig. [7] 
we compare the iso-contours corresponding to N iv/N m = 1 
(03.5 I/III or 03 V) for different wind-strengths, from thin (se- 
ries 'A') to dense winds (series 'F'). All models have LMC 
background abundances and [N] = 7.78. Since lower metallic - 
ities counteract the impact of wind-strength, corresponding iso- 
contours for SMC conditions vary to a lesser degree, whilst the 
variations are larger for Galactic conditions, in both cases by few 
hundreds of Kelvin. 

Generally, the impact of wind-density on the N iv/N ra emis- 
sion line ratio is quite strong, particularly for larger surface grav- 
ities. From Fig.|7]it is evident that for a fixed log g cooler r e ff are 
required to produce similar line ratios at higher wind densities. 
Basically, we can combine the reaction into three groups, namely 
A' to 'C models, 'D' models, and models with winds stronger 
than 'D'. The first group consists of low-M models with line 
ratios which are almost insensitive to wind-strength. Models be- 
longing to the 'D' series display a slight reaction, since the wind 
begins to affect the line emission (either from both lines or, for 
hotter stars, only from N iV/14058, see Fig. [3), and in a different 
way, so that the line ratio becomes modified. From the 'E' series 
on, we found larger differences, e.g., at log g = 4.0, of roughly 2 
and 4.5 kK for 'E' and 'F' models compared to low-M ones. 

In conclusion, objects with the same T e ff and \ogg but sub- 
stantially different wind-strengths are predicted to display signif- 
icantly different line ratios. Consequently, the assigned spectral 
types would be no longer monotonic in T e $ if there would be 
a large scatter in wind-strength. One has to note, however, that 
a certain T^Aogg pair also implies a certain luminosity class, 
and that the wind densities per luminosity class are only mildly 
varying. Thus, the monotonicity of a spectral type-r e ff relation 
might still be warranted for a specified luminosity class, as long 
as the corresponding luminosity class indicators (for the earliest 



NIV4058/NIII4640(strength): series A-F 



4.5 






X 


x 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X X 


1 

X 




X 


X 




X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


x jk) 




XXX 


X 


X 




X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


* 


X 






X 


XXX 


X 


X 


4.0 


— X 


* 


X 


X 


X 


X 


X 


X 


X 


X 


X 




*✓ 




X 


XXX 


X 


X 




X 


X 


X 


X 


X 


X 


X 


X 


X 


X 








X X 


X 


XXX 


X 


X 




X 


X 


* 


X 


X 


X 


X 


X 


X 








X 


X X 


X 


X X 








X 


X 


X 


X 


X 


X 


X 


X 








X 


X 


X 












X 


X 


X 


X 


X 


X 


X 
























3.5 


— x 


* 


X 


X 


X 




& 


% 


X 










log 




-14.00 


































































log Q 




-13.50 


































log Q 




-13.15 


































log Q 




-12.80 








X 


X 
























log 




-12.45 


































log 




-12.10 






3.0 






























1 









3.5x10* 4.0x10* 4.5x10* 5.0x10* 5.5x10* 

Teff (K) 

Fig. 7. Dependence of Niv/14058/Niil14640 on wind-strength 
(series A'-'F'), for LMC models with [N] = 7.78. The iso- 
contours correspond to a nitrogen emission line ratio of unity. 

O-stars, He n/14686) allow for a reliable classification (but see 
Sect. [8]). 

4.2.3. Impact of nitrogen abundance 

After investigating the impact of background metallicity and 
wind-strength, we now concentrate on the influence of nitrogen 
content. At first glance, this is somewhat surprising. Abundance 
should have a marginal or only small effect on line ratios, since 
it should cancel out as long as the lines are not too strong and the 
ionization equilibrium is not disturbed. This statement, however, 
is only true if the lines form in a 'typical' way, i.e., are simple 
absorption lines (e.g., He 1/1447 l/Heii/14541 used for the classi- 
fication of not too early O-stars). In the case considered here, 
however, the lines are formed by complex and different NLTE 
processes, and different abundances might have a different im- 
pact on the formation of both lines (see Sect. 13.3b . so that the 
ratio might become affected. Moreover, the variation of [N] in 
early type stars can be much larger than, e.g., the variation of 
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MW-NIV4058/NIII4640(strength): series E, log Q - -12.45 
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LMC-NIV4058/NIII4640(strengt.ri): series E, log Q = -12.45 




SMC-NIV4058/NIII4640(strength): series E, log Q 




4.5x10* 
Teff (K) 



Fig. 8. Dependence of Niv/l4058/Niii/t4640 on nitrogen abun- 
dance, for three different background metallicities: Galactic (up- 
per panel), LMC (middle panel), and SMC (lower panel). All 
models have log Q = -12.45 (series 'E'). The iso-contours cor- 
respond to a nitrogen emission line ratio of unity. 



Fee, which amplifies the effect. Thus, it is mandatory to test for 
the impact of nitrogen content on the N rv/N m emission line ra- 
tio. 

For this purpose, we display in Figj8]the response of the line 
ratio on different nitrogen abundances as present in our model- 
grid, for each Z (Galactic, LMC, and SMC background abun- 
dances). To allow for a fair comparison, all models belong to se- 
ries 'E'. Again, we display only contours with N iv/Nm — 1, as a 
representative value. Note that for each Z, the range of nitrogen 
content is different, from the corresponding baseline abundance 
to enrichments of roughly 1.3, 1.6, and 1.9 dex for MW, LMC, 
and SMC objects, respectively (see legend and Table[TJ. 



Table 3. Stellar and wind parameters of models used to compare 
synthetic nitrogen line profiles from fas twind and cmfgen. Th e 
grid is a subset of the grid presented by iLenorzer et al.l d2004). 
For details, see Paper I. Even entries provide those parameters 
from our fastwind grid models which reproduce best the H/He 
spectra from cmfgen. 



Model 


code 




R, 


l°g# 


log 2 






(K) 


(Re) 


(cgs) 




d2v 


CMFGEN 


46100 


11.4 


4.01 


-12.43 


E4740 


FASTWIND 


47000 




4.00 


-12.45 


d4v 


CMFGEN 


41010 


10.0 


4.01 


-12.75 


D4140 


FASTWIND 


41000 




4.00 


-12.80 


s2a 


CMFGEN 


44700 


19.6 


3.79 


-11.99 


F4540 


FASTWIND 


45000 




3.80 


-12.15 


s4a 


CMFGEN 


38700 


21.8 


3.57 


-12.15 


F3935 


FASTWIND 


39000 




3.50 


-12.15 



Notes. All models were calculated with v m \ c = 10 kms , and with 
four different nitrogen abundances, [N] = 7.78, 7.98, 8.38, 8.78. Wind 
strength parameter, Q, calculated in units of M yr~' , R Q and km s -1 . 



Generally, an enhancement of [N] at fixed log g shifts the iso- 
contours towards higher temperatures. We checked that this shift 
indeed is related to the different [N] -dependencies of the specific 
formation mechanisms of Niv/14058 and N 111/14640 rather than 
to a modified ionization equilibrium, which remains quite unaf- 
fected from variations in [N]. Moreover, the shift is quite similar 
for all Z, producing an increase of roughly 1 kK for a change of 
+0.2 dex in nitrogen abundance (at log g- 4.0; lower Ar e ff are 
found at lower log g). This is an interesting result because objects 
with considerable differences in [N] would display, for the same 
line ratio, i.e., spectral type, large differences in T e ff. If we, e.g., 
consider LMC models at log g = 4.0, a line ratio of unity (cor- 
responding to an 03 dwarf) would be obtained at temperatures 
differing by 6 kK if either [N] is at the baseline abundance (red, 
r eff * 42.5 kK) or if [N] = 8.18 (dark green, r eff ~ 48.5 kK), cor- 
responding to a typical enrichment of LMC stars (paper II)! Such 
a difference is much larger than the typical spread of temperature 
per spectral type, and we conclude that the present classification 
scheme for the earliest O-stars might be strongly biased by nitro- 
gen abundance. Only if the nitrogen content/enrichment would 
be rather similar for a given position in the r e ff-log g planfl this 
bias would not contaminate a spectral-type-r e jf relation, similar 
to our argumentation regarding the bias by wind-strength from 
above. 

5. Comparison with results from cmfgen 

In Paper I we compared our results for N m lines with those from 
the alternative code cmfgen, for a grid of models comprising 
dwarfs and supergiants in the T e g range bet ween 30 and 45 kK, 
and w ith 'old' solar abundances according to Grevess e & Sauvaj 
(Il998h . In the following, we examine the consistency also with 
respect to N iv and N v, concentrating on the hotter models (d4v , 
d2v, s4a, s2a, see Table [3}, but allowing for different nitro- 
gen abundances to check the predictions from the previous sec- 
tions. 

cmfgen models were computed with a modified photospheric 
structure following the approach from Santolava -Rev et a"D 



8 If, e.g., nitrogen would be enriched by rotational mixing alone and 
the initial rotational speeds were similar. 
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NIV4058 model series D, Galactic abundances NIV4058/NIII4640 model series E, Galactic abundances 
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Fig. 9. N iV/14058, switching from emission to absorption with Fig. 10. Equivalent width ratio N IV/14058/N ni/14640, as a func- 
increasing abundance. Model series 'D' (log<2 = -12.8), tion of nitrogen abundance. Model series 'E' (log Q — -12.45), 
logg = 4.0, compared to results from cmfgen, model d4v (at log g- 4.0, compared to results from cmfgen, model d2v (at T e s 



T eS « 41 kK, filled circles). 



: 46 kK, filled circles). 



dl997l) . smoothly connected to a beta velocity law. In our ap- 
proach, the Rosseland mean from the original formulation was 
replaced by the more appropriate flux-weighted mean. Several 
'exa ct' photospheric structures from tlusty 
showed exce llent agreement with our 



comparisons using exa 
(lHubenv&Lanz|[l99l 



method (see also Najarro et al.l 1201 lb . A turbulent velocity of 
10 kms~' was assumed when computing both the level popula- 
tions and line profiles. Our models account for the presence of 
H, He, C, N, O, Si, P, S, Fe, and Ni, totaling 3965 full levels 
(1319 super-levels) and * 70,000 lines. 

First, we examine the potential switch of NiV/14058 from 
emission to absorption when increasing the abundance. From 
Fig. 3] this effect should occur at lower wind-strengths and not 
too high r e ff, i.e., potentially for model d4v. Indeed, the pre- 
dicted effect is clearly present in the corresponding cmfgen mod- 
els, where ,44058 is in emission for [N] = 7.78 and 7.98, and in 
absoiption for [N] = 8.38 and 8.78 (filled circles in Fig.g). The 
total variation in equivalent width as a function of [N] compares 
very well to results from fastwind, which appear at somewhat 
higher T e s= 42.5 kK, as shown by the solid lines (based on our 
fastwind grid for log g- 4.0 and model series 'D'). Thus we con- 
clude that the predicted effect is more or less code independent. 

In Fig. [10] we compare the influence of the nitrogen abun- 
dance on the emission line ratio N iv/N m (calculated from EWs), 
by means of model d2v. Again, cmfgen (filled circles) predicts 
a similar effect as fastwind, i.e., the EW ratio increases with 
[N]. The actual values of these ratios, however, are quite differ- 
ent from those calculated by fastwind. Similar ratios are only 
reached at higher r e ff, because, as discussed below, fastwind 
predicts stronger N m and weaker N iv emission lines for hot ob- 
jects. Nevertheless, the '[N]-effect' is clearly visible in the cmf- 
gen models, not only for model d2 v but also for the others. Thus, 
the position of a certain emission line ratio in the r e ff-log g plane 
depends on the actual nitrogen abundance. As already discussed, 
this can lead to an ambiguity of the spectral-type-r e jf relation. 
We come back to this problem in Sect. Q 

In Appendix|A]we provide a detailed comparison of strategic 
H/He/N lines predicted by cmfgen and fastwind, for all models 
from Table [3] and for different [N]. fastwfnd spectra have been 
taken from our grid, to show that the provided resolution in p a- 
rameter space is sufficient in most cases. 



Let us first concentrate on the H/He spectra (upper panels 
of Figs. IA. ll IA.41 |A51 and lA.8b . which turn out to remain unaf- 
fected by typical variations in [N]. In most cases, there is an ex- 
cellent agreement between the H/He spectra from cmfgen and the 
closest or almost cl osest fastwind grid- model, even for the He i 
singlet lines (but see lNaiarro et al.l2006l) . The Stark-wings of the 
Balmer lines are well reproduced at a similar gravity. The largest 
discrepancies occur for model d2v where He lAAAAl 1 ,47 1 3 re- 
sult in a better fit at r e ff = 47 kK instead at the nominal value of 
46 kK, and for model s4a where log g — 3.5 produces better con- 
sistency than a model at log g — 3.6 which would lie closer to the 
nominal value of logg = 3.57. Finally, H a (fastwevd) shows less 
emission for model s2a, mostly because the closest grid model 
has a lower wind-strength of log<2 = -12.15 than the nominal 
value, log Q = -1 1.99. Figure lATl illustrates how the H/He spec- 
trum changes when is modified by +1000 K. In the consid- 
ered temperature range, the major reaction occurs in He i, pre- 
dominantly in the singlet lines (/L14387, 6678). 

The lower panels of Figs. lA~TllA.4llA.5l and lA.8l displav im- 
portant nitrogen lines for the four investigated models, at [N] = 
8.78 (ten times solar), whilst Figs. lAT2llA31IA.6l and lA"!71 displav 
these lines for [N] = 7.78 (solar). 

Nm lines have been already compared in Paper I, for 
[N] = 7.92 (and v m i c = 15kms _1 ), and these findings remain valid 
also at higher nitrogen abundance. 

For the dwarf models d4v and d2v, the emission lines at 
/1/14634-4640-4642 are stronger than predicted by cmfgen, and 
indeed stronger than displayed in Paper I, where the disagree- 
ment was less obvious. This discrepancy could be tracked down 
to a different value of v m ; c used in the present models. A lower 
value (10 vs. 15 km s _1 ) results in narrower resonance zones, and 
can lead to higher emission peaks0 Moreover, differences in the 
occupation numbers - when close to inversion in certain regions 
- become more pronounced, due to less averaging. For models 
s4a and s2a we note similar discrepancies, though at a tolerable 
level. 



9 A similar influence of v tur b has been seen, e.g., in syn thetic spectra 
of Br a when in emission due to photospheric processes (Naiar roet al.l 

EmTj). 
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As discussed in paper I, the Nm absorption line at /14097 
is slightly stronger in fastwind models at hot temperatures. 
Note that for d2v the profiles become identical when compar- 
ing with a model that reproduces the cmfgen Hei lines, i.e., for 
r eff = 47 kK instead of 46 kK. 

The quartet lines^ at AA4510 - 4514 - 4518 behave quite 
interestingly. Among the low [N] models, these lines are in ab- 
sorption only at d4v, whilst for the other models (higher T e g 
and/or higher M) they appear in emission. In a certain param- 
eter range (models d2v and s4a) and in analogy to NiV/i4058, 
these lines switch, for increasing [N], from emission to absorp- 
tion. Regarding these major trends, cmfgen and fastwind behave 
similarly, and the lines agree quite well for the dwarf models. 
For the supergiant models, on the other hand, fastwind predicts 
more emission/less absorption than cmfgen, in particular for s2 a. 
Thus, and due to their complex behavior, the quartet lines need 
to be treated with care when analyzing early O-stars. 

N iilI4003 and A4 195 (in the blue wing of Hen^4200 and 
additionally blended by Sim at hotter temperatures) show a 
mostly reasonable agreement. Nni/14195 matches almost per- 
fectly when predicted to be in emission, and [N] is low. The 
biggest discrepancy for /14003 is found in the low [N] model of 
s4a, where cmfgen predicts much more absorption. 

N iv and N v lines have not been compared so far. N iV/14058 
from fastwind shows always less emission than produced by cm- 
fgen, though for models s2a and for the low [N] model s4a the 
agreement becomes satisfactory. 

In combination with more emission from the N m triplet at 
hotter temperatures, this leads to a lower N iv/N m EW/line- 
strength ratio in fastwind (see Fig. [Toi l, i.e., effective tempera- 
tures deduced from this line-ratio alone would be higher com- 
pared to cmfgen results. This general trend is independent of 
abundance. 

Niv/13480 is in good agreement for most models, except 
for d4v where fastwfnd predicts more absorption. In contrast, 
NiV/16380 shows considerable differences, and for almost all 
models (except for the high [N] model of d4v) our grid predicts 
much more absorption. Since this line provided no difficulties 
when comparing with observations (neither in Paper II nor in 
the present investigation), we suggest that it might suffer from 
certain problems in cmfgen. We are currently working on this 
problem, which we think is connected to the desaturation of the 
N iv resonance line around 247 A at the base of the wind. 

Finally, the N v doublet /L14603-4619 compares very well in 
most cases, and only for the low [N] models of s2a and s4a we 
find more absorption by fastwind. 

Summarizing the major discrepancies, we note that in the 
early O-type domain fastwind produces more emission in the 
Nm triplet, less emission at NiV/14058, and mostly much more 
absorption at Niv/16380, compared to cmfgen. We now ask how 
such discrepancies might affect an abundance/r e ff analysis, and 
how one might proceed to diminish the impact of corresponding 
uncertainties. 

Since the H/He lines are in very good agreement, the stel- 
lar parameters should be constrained quite well, as long as He i 
remains visible. For the hottest stars with no or negligible He i, 
the situation might become more difficult, since r e ff needs to be 
constrained mostly from the nitrogen lines, and we would de- 
duce lower r e ff from cmfgen if concentrating on the N m triplet 
and Niv/14058 alone (see above). In turn, this would lead to 
rather different abundances. Fortunately, however, there is also 

10 Preferentially used by Mar tins et al.l ( l2012at) to derive nitrogen 
abundances in magnetic O-stars of late and intermediate spectral type. 



N v which is very sensitive on T e s (see below and Figs. IA.lllA.3l 
and I A. 4b . and seems to be quite code-independenO- In conclu- 
sion, r e ff for 'cooler' objects should be mostly determined from 
the He i/He n ionization balance. Potential discrepancies regard- 
ing N mXl4634 - 4640 - 4642 and N iV/14058 because of 'erro- 
neous' predictions can be easily identified then. For the hottest 
objects, on the other hand, Nv should obtain a strong weight 
when determining T e g. As long as either N iV/16380 or N iv/13480 
are observed, these lines might be used as abundance indicators, 
whilst any inconsistency of N iv/14058 will tell about the quality 
of that line. 

In Fig.QT]we now 'analyze' synthetic spectra from cmfgen 
by means of our fastwind grid, to check the overall consistency 
when accounting for all important lines. This is done for the high 
[N] model of d4v, which is a prototypical case, and by means 
of EW iso-contours in the r e jf-[N] plane (for the appropriate 
model series 'D' with log Q = -12.8, and a gravity of logg= 
4.0). Equivalent width for all analyzed lines have been mea- 
sured from the cmfgen spectra. These EWs are displayed as iso- 
contours with respect to our fastwind grid and should cross, in 
the ideal case of a perfect consistency between both codes, at one 
point corresponding to the model parameters, i.e., - 41 kK 
and [N] = 8.78. Because of the various discrepancies discussed 
above, this is not the case, but at least almost all iso-contours 
are close to each other in quite a confined region, marked by 
a black box. In particular, five out of eight lines cross around 
the point r eff = 41.5 kK / [N] = 8.55. At the original (cmfgen) 
values, r e (j = 41 kK and [N] = 8.78, on the other hand, there 
is only Nv/14603 and Nni/14003. Thus, a fastwind analysis of 
this model spectrum would yield lower nitrogen abundances, by 
roughly 0.2 dex (more on this below). 

Figure[TT]displays a number of interesting aspects. Most iso- 
contours have a parabola-like shape with a minimum at low [N] , 
where the 'left' branch with negative slope relates to increas- 
ing ionization fractions (the same EW at higher T e s implies a 
lower abundance), and the 'right' branch with positive slope to 
decreasing fractions. Because the right branch applies for lines 
from lower ionization stages (N m), and the left branch for lines 
from higher ones (N rv, N v), these different slopes allow to con- 
strain the actual parameter region quite well, even though both 
codes produce different line-strengths when compared at the 
same location in parameter space. 

The only feature which is quite outside the enclosed region is 
the N in triplet, which is significantly different in both codes (and 
stronger in fastwind, thus formed at lower T e g when compared 
to the cmfgen EWs). 

At the comparatively 'cool' temperature of model d4v, 
N V/L14603-46 19 becomes almost independent of abundance (al- 
most vertical left branch), and thus a very sensitive tempera- 
ture indicator. Indeed, it perfectly fits at the actual T e ff (which 
is also true for the He i/He n lines, see Fig. lA.U . Interestingly as 
well, the iso-contours for Niv/16380 and Niv/13480 occupy an 
almost identical region, i.e., give very similar results, though be- 
ing members of different spin systems (singlet and triplet system, 
respectively). Insofar, the information provided by both lines 
is similar, and the observation/analysis of either of these lines 
should be sufficient, which is fortunate since the 3500 A region 
is scarcely observed. 

Overall, when 'analyzing' [N] from the cmfgen spectrum and 
accounting for the rather fixed r e ff and log g values from H/He 
and N v, we would derive a lower or roughly similar nitrogen 

11 And not affected by X-rays from wind embedded shocks under typ- 
ical conditions, see Paper II. 
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Model series D, Galactic abundances 
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Fig. 11. Nitrogen abundance 'analysis' of cmfgen model d4v with [N] = 8.78. EW iso-contours for important nitrogen lines are 
displayed in the r e ff-[N] plane of our fastwind model grid (series 'D\ log Q = -12.8, logg= 4.0). The displayed iso-contours refer 
to EWs measured from the cmfgen spectra of model d4v. In the ideal case, i.e., if cmfgen and fastwind produce identical results, 
all iso-contours would cross at a single point located at T e ff- 41kK and [N] = 8.78. The black box marks the smallest region in 
parameter space where almost all lines (except for N ni/14640) are predicted at the measured value. In particular, five out of eight 
lines cross around the point T e s = 41.5 kK / [N] = 8.55, i.e., at a lower abundance (by roughly 0.2 dex) than present in the original 
cmfgen model. 



abundance from most lines, compared to the actual value used 
by cmfgen, except for the N m quartet lines which would imply a 
somewhat higher abundance. These results are also valid for the 
other models from Table[3] We thus conclude that in those cases 
where T e g can be additionally constrained, fastwind analyses of 
early O-type stars will yield mostly lower nitrogen abundances 
than analyses performed by means of cmfgenFI 



LMC sample analyzed in Paper II contained only a few early- 
type objects (BI237, BI253, Nl 1-026, Nll-031, and Nl 1-060). 
These stars, biased towards 02, did not cover all spectral types 
and luminosity classes within the classification scheme. To allow 
for a larger sample, we adde d the earliest s t ars from the analy- 
ses of LMC/SMC objects by iMassev et al.1 (120041 12005L I2009L 
hereafter Mas04, Mas05, and Mas09). 



6. Analysis of a sample of LMC/SMC early-type 
O-stars 

So far, we mainly provided theor etical predictions on t he N rv/ 
N in emission line ratio as used by Wal born et al.l (120021) to clas- 
sify the earliest O-stars, and concentrated on the impact of vari- 
ous parameters. For testing our predictions, we need to confront 
these results with the analysis of a significant stellar sample. The 

12 At the present stage of knowledge, we do not know which code 
does a 'better' job in reproducing reality. 



6. 1 . The stellar sample 

In particular, we selected objects with spectral types earlier than 
05. Table|4]gives a complete list of all early O-type stars consid- 
ered in the following, along with their spectral type and galaxy 
membership. The final sample consists of seventeen stars, four- 
teen from the LMC and three from the SMC. Thirteen objects 
have been drawn from Mas04 and Mas05, plus one object from 
Mas09 (the only early-type star of that sample). Two of the ob- 
jects analyzed in Paper II (BI237 and BI253) were also studied 
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Table 4. Sample stars used within this study, along with galaxy 
membership and spectral type. 



Star 


Cross-IDs 


Galaxy 


Spectral Type 


AV 177 


- 


SMC 


04 V((f)) fl 


AV435 


- 


SMC 


03 V((f ))" 


NGC 346-355 


NGC 346 W3 


SMC 


ON2 111(f)* 


LH81:W28-5 


- 


LMC 


04 V((f+)) fl 


LH 81:W28-23 


- 


LMC 


03.5 V((f+)f 


LH 90: ST 2-22 




LMC 


03.5 III(f+) fl 


LH 101:W3-24 


ST 5-27 


LMC 


03.5 V((f+)) fl 


y u in! -WTI 1 Q 
LI 1 1U1 . W j-iy 


Ol J-jl 




UZ 11 


R136-018 




LMC 


03 111(f)" 


R136-040 


R136a-535 


LMC 


02-3.5 V" 


Sk-65° 47 


LH 43-18 


LMC 


04If 


Sk-67° 22 


BAT 99- 12 


LMC 


02 If /WN5 C 


BI237 




LMC 


02 V((f ))" 


BI253 




LMC 


02 V((f )) d 


Nl 1-026 




LMC 


02 III(f ) d 


Nll-031 


P3061/LH 10-3061 


LMC 


0N2 III(f f 


Nl 1-060 


P3058/LH10-3058 


LMC 


03 V((f )Y 



(a) Mas05; {b) IWalborn et al.l ll2004h : (c) ICrowther & Walbornl (1201 lh : 
(rf) lMokiemetalJ ( l2007iJ) and references therein 

Notes. The horizontal line separates stars analyzed within 
this work (drawn from Mas04/05/09) from stars previously 
analy zed in Paper II. I d entific ations are as follows: "AV" 
from lAzzopardi & Vi^nea ul (Il982h; " BAT" from Brevsache r et al J 
dl999h : "BI"from iBrunet et al.l ( I1975T): "LH" from iLuckd d 19721) 
excep t "LH 10-305 8" that is from lwalborn et alj d2002L [20041): "Nl 1 
from 
from 



Evans et al.l J200i 
Parker et al} (fl99: 



NGC" from Massev et al. (1989); "P" 
_ R136" from iMassev & Hunted il998h: 
R136a" from Malu muth & Head <1994l) : "Sk" from ISand uleak 



by Mas05. For these objects, as well as for the remaining sam- 
ple members (N 1 1 -026, N 1 1 -03 1 , and N 1 1 -060), we refer to our 
analysis from Paper II. 

From the original subsample of early-type stars studied by 
Mas04 and Mas05, we selected representative objects. Owing to 
various reasons, eight stars have been discarded, (i) Three dwarfs 
(R136-033, R136-040 and R136-055) did not show any trace of 
either N ni/U4634 - 4640 - 4642 or N iV/14058 in their spectra, 
and Massey and co-work ers could not classify them according to 
the IWalborn et all {2002) scheme. Therefore, we analyzed only 
one prototypical example for these stars, R 136-040. (ii) From the 
four early giants compiled in Mas05, we discarded two of them, 
R136-047 (02 111(f)) and LH 64-16 (02N 111(f)). The FOS 
data used for the R136 stars, see Sect. 16.21 suff ered from an inter- 
mitte nt behavior of some of the FOS diodes ( Masse v" & Hunterl 
IT998h . which could result in the appearance of spurious fea- 
tures contaminating the spectra. Owing to a rather bad quality 
of the spectra, we were not able to obtain a satisfactory fit for 
R136-047, and we decided to discard it from the analysis. This 
problem was not met for the remaining R136 stars used within 
this work (R136-018 and R136-040). The other giant discarded, 
LH 64-16, shows similar severe discrepancies as we had found 
in Paper II for another star of this class, Nl 1-031. Owing to ex- 
treme difficulties to fit either N m/N iv or N iv/N v lines, we ex- 
cluded it from the present analysis, (iii) Regarding supergiants, 
we selected three out of seven stars. For typical conditions (M 
> lO^Moyr 1 , v M ~ 3000 kms" 1 , and R t ~ 15 R G ), the corre- 
sponding wind-strength, log Q > -12, is well outside of our 
model-grid (Table[TJ, and we concentrated on representative ob- 



jects since analyzing all of them would have been extremely 
time-consuming, even with a fast code such as fastwind. 

In summary, our sample comprises 9 dwarfs, 5 giants, and 3 
supergiants. The source of the corresponding spectral types for 
each object is listed in Table |4] 

6.2. Observational data 

Spectra for the 12 early O-stars from Massey et al. were taken 
in three ranges, the UV, the blue and the red visual band. Only 
optical data were used within this work. 

(i) For the bulk of the stars, these data were obtained at the 
Cerro Tololo Inter-American Observatory (CTIO) 4m Blanco 
telescope with the RC spectrograph, during January 1999 (P.I. 
P. Massey). The blue and red spectra cover the ranges 3750- 
4900 A and 5400-7800 A, respectively. The spectral resolution 
was 1 .4 A for the blue and 2.8 A for the red band, with a S/N be- 
tween 200 to 560 per resolution element, corresponding to about 
4 pixels, in the blue (measured at 4500 A), and a typical S/N of 
about 150 per resolution element in the red. 

For LH 10LW3-24, additional spectra for the H a region 
were taken by means of the Hubble Space Telescope (HST) us- 
ing STIS/CCD (as part of programme GO-9412, P.I. P. Massey), 
to amend the ground-based observations which were heavily 
contaminated by strong nebulosity. These data cover 6300- 
6850 A, with a spectral resolution of 0.84 A and S/N of 50. 

(ii) For the R136 stars, we used two different data sets ob- 
tained with the HST One was taken with the G400 FOS/RD 
configuration under GO-6417 (P.I. P. Massey), with a resulting 
wavelength coverage from 3250 to 4820 A, 3 A resolution and 
a S/N of 60 at 4400 A. Additional observations were taken with 
STIS/CCD under GO-7739 (P.I. P. Massey). Blue observations 
covered the wavelength range 4310-4590 A, at a S/N of 100 for 
a spectral resolution of 0.4 A. Observations centered around H a 
were made with the same setup as for LH 101:W3-24 described 
above. Since the STIS spectra have better S/N and resolution, 
we used them preferentially for our analysis. Because of their 
limited wavelength coverage, excluding the spectral region that 
comprises the most important nitrogen lines, we nevertheless 
had to consider the FOS data, see Figs. IB.ll and lB.5l 

(iii) Spectra for NGC 346-355 were collected with the Clay 
6.5m (Magellan II) telescope at the Las Campanas Observatory, 
using the Boiler & Chivens Spectrograph (P.I. P. Massey). The 
coverage was 3410-5040 A with a S/N of 600 at 4500 A in the 
blue, and from 53 15 to 6950 A in the red, achieving a S/N of 400 
at 6500 A. For both ranges, the spectral resolution was 2.4 A. 

All spectra were reduced with irafQ and we performed ad- 
ditional re-normalizations for different wavelength ranges within 
this work. 

As outlined in Sect. 16.11 our present sample includes two 
objects (BI237 and BI253) which have been studied both in 
Paper II and by Mas05. The UVES spectra used in the former 
analysis have a considerably higher resolving power, 40,000, 
compared to a resolving power of 3,000 as provided by the RC 
spectrograph used by Mas05. After comparing both data sets 
amongst each other and with our fits from Paper II, it turned out 
that the nitrogen lines are quite similar, whilst the cores of the 
hydrogen (except for H„) and He n lines are somewhat deeper in 



13 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 



Research in Astronomy, Inc.. 
National Science Foundation. 



under cooperative agreement with the 
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the UVES data, more than expected from the higher resolution. 
Because of the high quality of the UVES data, we kept the pa- 
rameters as derived in Paper II, though a re-analysis based on the 
Mas05 data would have provided rather similar values. 

6.3. Analysis 

Stellar and wind parameters were derived following the method- 
ology outlined in Paper II, and are listed in Table Comments 
on individual objects and corresponding H/He/N line fits are pro- 
vided in Appendix IB1 

In brief, we proceeded as follows. In accordance with 
Paper II and the investigations by Massey and collaborators, 
we assumed unclumped mass-loss throughout the analvsisPl 
First, we determined v sin i for each object using the Fourier 
method (lGravl[T976 )F^ Subsequently, we roughly constrained 
the stellar/wind parameters as well as the helium and nitrogen 
abundances by means of our model-grids for LMC and SMC 
background metallicities. For deriving T e ff, we used both the he- 
lium and nitrogen ionization balance when possible, i.e., when 
nitrogen lines from at least two ionization stages were visible For 
the cooler objects of our sample (T e ff < 44 kK), we mainly relied 
on the helium ionization balance, and used nitrogen as a consis- 
tency check. For hotter objects, where He n/14541 begins to loose 
its sensitivity to T e ff, we gave larger weight to the nitrogen bal- 
ance, as long as a reasonable fit to the He i/He n lines could be 
maintained. Note that any r e ff/[N] degeneracy (as present, e.g., 
for the line-strength ratio of N iv/N m) is broken by the absolute 
strengths of the lines, and the fact that we have usually more 
than two lines at our disposal, which react quite differently on 
changes in r e ff and [N] (cf. Fig. fTTb. Thus, as long as lines from 
at least two nitrogen ionization stages are available, T e ff and [N] 
can be determined in parallel from fitting these lines (as long 
as the other stellar/wind parameters can be inferred from inde- 
pendent diagnostics). If nitrogen lines were not available (too 
weak), we either relied on He i/Hen alone (LH 101:W3-24 and 
AV 177), or, when even Hei was no longer visible, we settled 
for a lower limit on T e g (R136-040). We refer to AppendixlBlfor 
a detailed discussion on the specific diagnostics used (see also 
Table[5]), and on particular problems for individual stars. 

Gravities were inferred from the wings of the Balmer lines, 
and log Q from fitting H a and He H/14686, with nitrogen lines 
serving as a consistency check for both quantities. Because of 
the rather low resolution and the mostly high v sin i values, we 
were able to constrain the 'macro-turbulence', v mac , for only two 
objects (LH 81:W28-5 and LH 90:ST 2-22). For the remain- 
ing ones, our fits were acceptable without any need for extra- 
broadening. 

At this point, we adopted terminal velocities as derived by 
Massey et al. from UV lines, updating M to preserve Q if nec- 
essary. For the bulk of the stars, we kept the velocity field ex- 
ponent used within our model-grid, /3 = 0.8, because of rea- 
sons of consistency, since also Massey et al. used this value. 
Only for LH 8LW28-23 we adopted /3 = 1.0, to better repro- 
duce He H/14686 (see Fig. IB. 2b . For those objects with Hq, and 
He ri/14686 in emission (LH 10LW3-19, Sk-67° 22, and Sk-65° 
47), it would have been also possible to derive /3 from the pro- 
file shape. Because of the quite good fit-quality for both lines 



14 Regarding the impact of clumping on the synthetic nitrogen spectra, 
see Paper II. 

15 As implemented and tested in the OB-star range by 
ISimon-Dfaz et al. (2006h and lSimon-Dfaz & Herrerol (l2007h . 



already with = 0.8 (see Figs. |B77] and US, we had no 
real reason to change this value though. 

Final parameters were derived from a grid of much higher 
resolution around the initial constraints. Stellar radii were calcu- 
lated from My as provided by Mas04/05/09 and synthetic fluxes, 
with a corresponding final update of M. Following our experi- 
ence for the earliest O-stars from Paper II, we used a micro- 
turbulent velocity, v m ; c = 10 kms~', for all objects considered 
here. 

Errors on stellar/wind parameters are adopted following 
Paper II. In brief, we estimate an uncertainty of +1 kK in 
Teff (for those objects with a unique solution, see Table |5), 
+0.1 dex in logg, and +0.05 dex in log /?,, adopting a typica l 
error of +0.25 mag in M v (see Eq. 8 in Repo lust et all 12004). 
Together with the errors for r e jf, this adds up to an uncertainty 
of +0.11 dex for log L/L Q . Larger errors, plus/minus a factor 
of two, are present for M, when H a is in absorption[3 For ob- 
jects with mass-loss indicators in emission, the error is some- 
what lower. Errors on v«, are on the order of +100 kms -1 , taken 
from Massey et al., and we estimate the errors on v sin i and 
Vmac as +10 to 20 km s _1 . A rough estimate on the error of Yu e is 
+0.01 to 0.02. Regarding [N], we adopt a conservative value of 
+0.15 to 0.20 dex to account for the dependence on stellar and 
wind parameters. 

6.4. Comparison with results from Massey et al. 

In the following, we briefly discuss overall differences between 
the stellar/wind parameters derived within this work and by 
Massey et al. for overlapping objects, also regarding BI237 and 
BI253 which were already analyzed in Paper II. For details, see 
AppendixlBl 

For the majority of objects, where we still could use the 
He i/He n balance as a primary temperature indicator (but also 
accounting for nitrogen, see Sect. 16.3b . we derived slightly cooler 
r e ff, on the order of 0.5 to 1 kK. This is not too disturbing, how- 
ever, as one often finds systematic differences of this order when 
different elements are used to derive T e s for hot stars, as, e.g., 
in the case of using Si vs. He for early B-/late O-type stars (e.g., 
iHunter et alj |2007). The maximum difference amounts to AT e f{ 
= -2 kK for LH 8 1 :W28-23, on the margin of the adopted errors. 

For a few objects (BI237, BI253, NGC 346-355, R136-018, 
and Sk-67° 22) where we needed to exploit the nitrogen ion- 
ization balance to break the helium degeneracy, we inferred 
considerably hotter temperatures (Ar e ff ~ 2-7 kK). Admittedly, 
and owing to the restricted quality of the data for some of 
the earliest objects with weak Nm triplet emission and very 
weak He 1/1447 1 , the effective temperatures derived might by af- 
fected by uncertainties due to noise and continuum placement. 
Nevertheless, the fact that also here the differences are system- 
atic suggests that the solution to this problem may not be purely 
observational. 

We also derived quite similar values for log g, except for the 
five objects explicitly mentioned above, which indicated higher 
gravities (mostly because of the higher r e ff), with a maximum 
difference of roughly 0.25 dex for BI253. Our M values are sys- 
tematically lower, mostly because of the lower T e g. For those 
stars with increased T^, the reduction is caused by lower R* 
and/or higher /3, where the latter effect is particularly strong 
for BI253 which displays the largest difference, AlogM = 
-0.5 1 dex. Helium abundances agree quite well except for those 



16 In this case, one is unable to derive [} from li ne-fitting, and the un - 
certainty in /3 reflects in quite a large error for M (Ma rkova et alj2004h . 



14 



J.G. Rivero Gonzalez et al.: The earliest O-stars 
Table 5. Fundamental parameters for the early O-star sample, assuming unclumped mass-loss. 



Star 


Spectral type Nitrogen 






log gtrue 




logL, 




log 2 p 


Voo 


v sin 


(' v mac 


*He [N] 






T eff diag. 


(kK) 


(cgs) 


(cgs) 


(Ra) 


(Le) 






(km s 


) (kms~ 


') (kms- 


') 


Dwarfs 


BI253" 


02 V((f )) 


Niv/Nv 


54.8 


4.18 


4.20 


10.7 


5.97 


1.53 


-12.61 1.21 


3180 


230 




0.08 7.90 


BI237" 


02 V((f )) 


Niv/Nv 


53.2 


4.11 


4.12 


9.7 


5.83 


0.62 


-12.98 1.26 


3400 


140 




0.09 7.38 


R136-040 


02-3.5 V 




>51.0 


4.01 


4.02 


10.3 


5.81 


1.93 


-12.53 0.80 


3400 


120 




0.08 6.90 


Nl 1-060" 


03 V((f )) 


Nm/Niv/Nv 


48.0 


3.97 


3.97 


9.5 


5.63 


0.51 


-12.92 1.26 


2740 


68 


40 


0.12 8.20 






Niv/Nv 


51.0 


4.10 


4.10 


9.2 


5.71 


0.48 


-12.92 








8.15 


AV 435* 


03 V((f )) 


N in/N iv 


46.0 


3.90 


3.91 


13.8 


5.88 


0.21 


-13.15 0.80 


1500 


110 




0.10 7.58 


LH81:W28-23 


03.5 V((f+)) Nm/Niv/Nv 


47.0 


3.80 


3.82 


10.0 


5.65 


1.57 


-12.53 1.00 


3050 


146 




0.25 8.40 


LH 101:W3-24 03.5 V((f+)) 




47.0 


4.00 


4.01 


8.1 


5.46 


0.27 


-13.00 0.80 


2400 


120 




0.10 7.78 


LH81:W28-5 


04 V((f+)) 


Nm/Niv/Nv 


44.0 


3.80 


3.81 


9.8 


5.51 


1.08 


-12.60 0.80 


2700 


120 


80 


0.15 8.38 


AV 177* 


04V((f)) 




44.0 


3.80 


3.85 


8.8 


5.42 


0.23 


-13.19 0.80 


2650 


220 




0.15 7.78 


Giants 


Nl 1-031" 


ON2 III(f ) 


N ni/N iv 


47.8 


3.95 


3.95 


13.4 


5.92 


2.02 


-12.64 1.08 


3200 


71 


60 


0.11 7.83 






Nrv/Nv 


56.0 


4.00 


4.00 


12.2 


6.12 


2.20 


-12.54 








8.30 


NGC 346-355* 


ON2 IIKf) 


Nrv/Nv 


55.0 


4.00 


4.01 


12.0 


6.08 


2.50 


-12.39 0.80 


2800 


140 




0.10 8.10 






N in/N iv 


51.0 


4.00 




12.5 


5.98 


2.00 


-12.51 








7.98 


Nl 1-026" 


02 111(f) 


N m/N iv/N v 


49.0 


4.00 


4.00 


11.3 


5.82 


1.56 


-12.63 1.08 


3120 


72 


60 


0.10 7.80 






Niv/Nv 


52.0 


4.10 


4.10 


11.0 


5.89 


1.49 


-12.63 








7.75 


R136-018 


03 ni(f* ) 


N m/N iv/N v 


47.0 


3.90 


3.92 


14.2 


5.95 


1.68 


-12.76 0.80 


3200 


180 




0.10 8.18 


LH 90: ST 2-22 


03.5 HI(f+) 


N m/N iv/N v 


44.0 


3.70 


3.71 


18.8 


6.08 


4.56 


-12.36 0.80 


2560 


120 


80 


0.15 8.58 


Supergiants 


Sk-67° 22 


02 If /WN5 


Nm/Niv/Nv 


46.0 


3.70 


3.74 


12.4 


5.80 


15.00 


-11.60 0.80 


2650 


200 




0.30 8.78 


LH 10LW3-19 


02 If 


Nm/Niv/Nv 


44.0 


3.90 


3.91 


25.3 


6.34 


20.86 


-11.97 0.80 


2850 


180 




0.10 8.18 


Sk-65° 47 


04If 


N m/N iv 


40.5 


3.60 


3.62 


19.8 


5.98 


11.00 


-11.89 0.80 


2100 


160 




0.12 7.78 



'"' Stellar/wind parameters derived in Paper II; (b) SMC star; w in units of 10 6 M Q yr 1 

Notes. logg tmc is the surface gravity corrected for centrifugal effects. For each star, the additional (or even primary) r c ff diagnostics used in 
parallel with the He i/He n ionization equilibrium is indicated. For four stars, two parameter sets are provided, because we had problems to reach 
a simultaneous fit for all considered lines from three nitrogen ionization stages. For Nl 1-026 and Nl 1-060 the first entry is our preferred one, and 
the second entry is a (disfavored) alternative, whereas for Nll-031 (see Paper II) and NGC 346-355 we consider both solutions as possible. For 
errors, see text. 



few stars with a larger change in Tgg, where we had to adapt 
Ym to preserve the fit quality of the helium lines. The largest 
difference found amounts to about 0.05 in Fee for LH 90:ST 2- 
22. Finally, most of our v sin i are in good agreement with the 
Massey et al. values, and only three objects displayed substantial 
differences, with the largest one of about 30 km s~'for BI253. 

One of the major implications of our re-analysis of the ear- 
liest O-stars as considered by Massey et al. can be stated al- 
ready now: We do not find the pronounced degeneracy of the 
N iv/N in emission line ratio as claimed by Mas05 (for details, 
see Sect. 17.3b . For instance, these authors questioned the mono- 
tonicity of the classification scheme with respect to T e g , because, 
among other problems, they inferred similar T e ff for two dwarfs, 
LH 10LW3-24 (03.5) and BI237 (02), which on the other hand 
displayed considerable differences in their Niv/Nm emission 
line ratios. These r e ff were derived by a pure H/He analysis. By 
exploiting the nitrogen ionization balance, we are now able to 
break this degeneracy, since, e.g., we find a considerably hotter 
temperature for BI237 (see table |5), similar to the temperature 
derived for BI253, the other 02 dwarf from the sample. 

7. Discussion 

7. 1 . Nitrogen abundances 

Though not the major topic of our present work, let us briefly 
comment on the nitrogen abundances derived within this analy- 
sis. Globally, these are consistent with our results from Paper II. 



The bulk of the LMC stars displays a considerable enrichment, 
stronger than expe cted from tailored evolutionary calculations 
dBrott et al.ll201 lb), and supporting the idea of an efficient mix- 
ing during very early phases of O-star evolution. We confirm the 
tight correlation between the derived helium and nitrogen con- 
tent. 

Only three dwarfs (note that only one prototypical case of 
nitrogen-weak objects, Rl 36-040, has been analyzed) seem to 
be unenriched, visible already in the absence of nitrogen lines in 
their spectra. The star with the highest enrichment within the 
sample is Sk-67° 22 (02 IF/WN5), and the derived [N] lies 
roughly 1 .9 dex above the LMC baseline abundance, larger than 
any of the values found in Paper II. In analogy, also its helium 
content is extreme, and both findings support the evolved nature 
indicated by its 'slash' designation. 

Although the nitrogen content (in absolute numbers) of the 
few SMC sample members lies below that of most of our LMC 
stars, all of them are strongly enriched, by more than one dex 
above the SMC nitrogen baseline abundanceQ The largest en- 
richment (by -1.6 dex, [N] = 7.98. . . 8.10) was found for NGC 
346-355, in agreement with its 'ON' designation. A very simi- 
lar abundance, [N] = 7.92, together with si milar stellar param- 
eter s, has been previo usly derived by both Bo uret et al.l {2003) 
and lHeap et al.l d2006l) for this star, see Appendix IBl 



17 According to lHunter et al](l2007h . [NWiine = 6.5 
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Spectral Type 

Fig. 12. Effective temperatures as a function of spectral type, for 
LMC and SMC O-stars analyzed within this work and Paper II 
(Nl 1-031 and R136-040 discarded, see text). For NGC 346-355 
(ON2 III(f ), SMC), T e g according to our 'average' solution, see 
text. Black symbols: LMC objects; blue symbols: SMC objects; 
green symbol: Sk-67° 22 (02 If /WN5, LMC). Squares refer 
to supergiants, triangles to (bright) giants, and circles to dwarfs. 
Filled symbols correspond to objects with [N] > 8.0, and open 
ones to objects with [N] < 8.0. The dashed-dotted (red) line dis- 
plays a least square linear and quadrati c fit to all dwarfs, ac cord- 
ing to Eq.Q] The T e g calibration from Martins et al. (2005]) for 
Galactic O-dwarfs (dashed) and supergiants (dotted) is displayed 
for comparison. For clarity, some objects have been slightly 
shifted horizontally. 



7.2. Effective temperatures vs. spectral types 

In Fig. Q~2] we display the derived effective temperatures as a 
function of spectral type for our present stellar sample (Table|4|, 
augmented by all (LMC-) stars later than 04 from Paper II to 
extend the sampling toward cooler types analyzed in a uniform 
way. 

Because it was not possible to assign a specific spectral type 
to Rl 36-040 (see above), this star is not contained in Fig. [12] We 
also discarded Nl 1-031 (ON2 III(f")) from this diagram, owing 
to severe problems in the determination of its T e g. In Paper II 
we were not able to derive a unique effective temperature from 
He i and all three nitrogen ionization stages in parallel, whilst 
using either Hei/Nm/Niv or Niv/Nv resulted in a difference 
of Ar e ff= 8 kK, which is too large to allow for further con- 
clusions. Nevertheless, our sample contains another ON2 giant, 
NGC 346-355 from the SMC. Also for this object, we were only 
able to obtain a hotter (based on N iv/N v) and a cooler (based 
on N m/N iv) solution, but the difference is much smaller, Ar e ff = 
4 kK. In Fig. [12] we assigned a mean value, T e g= 53 kK (ac- 
counting for larger errors than typical, on the order of 3 kK) , to 
remain unbiased from a somewhat subjective view, but note that 
both the hot and the mean temperature are considerably higher 
than what we would have derived using just the He n/He i lines 
alone. Further comments on the ON2 stars are given at the end 
of this section. 

To check the impact of [N] on the T e ff estimates for a given 
spectral type, in Fig. [12] we denote objects with nitrogen abun- 
dances above and below [N] = 8.0 by filled and open symbols, 
respectively. This threshold has been chosen according to our 
findings from Paper II, roughly separating LMC objects with 



mild or typical enrichment from those with unexpectedly strong 
enrichment. We note that this value would be too high for SMC 
objects if one would be interested in displaying the actual enrich- 
mentQbut for our purpose of testing and comparing the impact 
of [N], only the absolute value and not the enrichment is decisive 
(see Sect. 14331 . 

The complete LMC sample considered in Fig. Q~2]covers 26 
stars spread over (almost) the full range of spectral subtypes, 
comprising 16 dwarfs, 7 (bright) giants, and 3 supergiants. On 
the other hand, only three early SMC objects (two dwarfs and 
one giant) have been analyzed, preventing firm conclusions. 

Before we concentrate on the earliest types, we highlight 
some general trends. From a first glance, and even though the 
number of supergiants and (bright) giants within our sample 
is lower than the number of dwarfs, it is obvious that T e ff in- 
creases from supergiants to d warfs for all spectral types, similar 
to Ga lactic conditions (e.g., Rep olust et al.ll2.Q04l iMartins et alj 
2005). At least for types later than 03.5, giants and supergiants 
are about 1 kK and 4 kK cooler than dwarfs, respectively. The 
latter difference is simi lar to what has been f ound in previous 
studies, e.g., by Mas05, Mo kiem etall (l2007ah . and Mas09. For 
the earliest types, on the other hand, this difference becomes 
much larger, about 10 kK at 02. 

We can also clearly distinguish different behaviors of the 
spectral-type-r e ff relation. For luminosity classes V and II-III, 
later types follow a linear trend, whilst the increase in T e g is 
much steeper for the earliest ones. On the other hand, the few 
early supergiants of our sample seem to follow a linear trend 
with a slope similar to cooler objects. 

At least for the dwarfs, we are able to provide a typical rela- 
tion, when ignoring any differences in Z and [N]. Using a linear 
and a quadratic least-square fit for the objects later/including and 
earlier than 04, respectively, we find 



T e ff - 



51.64- 
70.87 ■ 



1.94 x ST 
10.29 xST + 0:, 



if ST > 4 
XST 2 if2<ST<4, 



(1) 



where ST is the spectral type for O-dwarfs, and T e g is expressed 
in kK. At the present state of knowledge, this relation might be 
applied to LMC stars only, since the low number of analyzed 
SMC-dwarfs as well as other arguments (see Sect. 17. 51 prohibit 
an application for stars from the SMC. 

For comparison, we show in Fig{12]the observed spectral- 
type-r e jf relati on for Galact i c dwa rfs (dashed) and supergiants 
(dotted) from Ma rtins et al.1 (120051) . For the dwarfs, there is a 
typical offset of roughly 1 kK, whereas the (cooler) 02 and 04- 
supergiant seem to follow the Galactic calibration. 

A similar comparison betwe en their LMC sampl e and 
Galactic dwarfs was provided by iMokiem et al.l (12007 al) . who 
found a somewhat larger offset by ~ 2 kK. This difference is 
caused by lower T e s as resulting from the updated version of 
fastwind for later spectral types (see Paper II). 

Reassuringly, the scatter in the spectral-type- T e g relation for 
objects later than 04 is small (~1 kK), since the spectral types 
(as well as our primary T e g indicator) rely on the He i/He n line 
strength (or E.W.) ratio, which is a rather monotonic function of 
T e ff for a given luminosity class (controlling gravity and wind- 
strength) and background abundance{3 The only outlier is Nl 1- 



18 Indeed, all SMC objects in Fig. [T2] appear with open symbols, de- 
spite their strong enrichment. Because of the lower baseline abundance, 
a threshold of [N] = 8.0 corresponds to an enrichment of 1.5 dex, which 
is very unlikely to occur except for extreme objects such as NGC 346- 
355 located close to this value. 

19 All our cool objects are LMC stars. 
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065, an 06.5 dwarf. The analysis of this star was quite diffi- 
cult, and for the determination of T e g (and other parameters) 
we used also the N m/N iv ionization balance, in parallel with 
He i/He n. Nevertheless, we found considerable problems in fit- 
ting the emission at N m/M4634 - 4640 - 4642 together with the 
pronounced absorption at N iv/14058 and the remaining nitrogen 
lines. Even though this problem became somewhat relaxed by 
allowing for wind-clumping, it seems that there was indeed a 
problem in the analysis of this star, and that it might be cooler. 

Our sample includes two 04 dwarfs, one from the LMC 
(LH 8 1 : W28-5) and the other from the SMC (AV 177). For these 
stars we derive similar T e ff and log g, see Table [5] though the 
SMC star ( because of lower Z and weaker winds) should be 
hotter (e.g., iBouret et al.ll2003l iHeap et al.ll2006l iMokiem et aT] 
2006). Mas04/05, who came to the same general conclusion, 
found from a pure H/He analysis the same T e g for AVI 77 
(SMC), whilst they derived an even higher T eff for LH 8LW28- 
5 (LMC). We note that the present spectrum of the latter shows 
weak emission in NiV/14058 (Fig. IB.4I ). so that the star might 
be alternatively classified as 03.5. Since also the gravities are 
quite low for lc V stars, the situation for both objects remains 
somewhat unclear, even though the derived parameters are fully 
consistent with the observed N iv and N m lines. 

We concentrate now on the impact of different parameters 
such as [N], Z, and M, on the spectral-type-T^ff relation for the 
earliest O-types (02-03.5), for which the spectral classification 
depends crucially on the N iv/N m emission lines, whilst our T e ff 
diagnostics includes additional lines from those ions as well as 
from N v when visible. To allow for an easy understanding of the 
following, we summarize the results from our parameter study 
in Sect. [4] For a given emission line ratio (i.e., spectral type), 
the derived T e ff should increase with [N] and log g, and decrease 
with Z (at least for higher T e s, see Fig. [6j and M (more precisely, 

log 6). 

Thus, the general r e ff-difference between dwarfs, (bright) gi- 
ants, and supergiants can be easily attributed to differences in 
log g and log Q, which increase and decrease for increasing lu- 
minosity class, respectively. In this way, both effects add up, 
leading to r e ff(lc V) > r e jf(lc III) > TWGc I), at least for com- 
parable [N] and Z. Regarding this general trend, a classification 
in terms of nitrogen (early O-types) and helium (late O-types) 
provides similar effects. 

• For the two LMC 03.5 dwarfs (LH 8LW28-23 and 
LH 10LW3-24), we find the same r eff , though logg dif- 
fers by 0.2 dex. The effect produced by the larger [N] of 
LH 8LW28-23 (A[N] * 0.6 dex compared to LH 101 :W3- 
24, implying a shift towards higher T e g) is counteracted by 
both its denser wind and a lower logg. Anyway, it is not 
clear whether LH 8LW28-23 is correctly designated as a 
dwarf. Indications of a giant nature are its low surface gravity 
(logg = 3.8), the wind-strength, and the trace of a P-Cygni 
profile at He H/14686 (Fig |B.2t . If this would be the case, the 
inferred T e ff might be too high for this luminosity class when 
compared to the 'other' 03.5 giant, LH 90:ST 2-22. 

• The two 03 dwarfs, Nl 1-060 from the LMC and AV 435 
from the SMC, show different T eff . Astonishingly, the LMC 
star is hotter than the SMC one, by 3 kK, contrary to what 
might be expected. Here, the [Nl effect outweighs the corre- 
sponding one associated to 2P1 whilst differences in \ogg 



20 In cooler stars, the He i/He n ionization balance used for classifica- 
tion is almost only affected by backgro und abundance (less Z, less EUV 
line-blocking), see Re polust et al.l fe004l) . 



and wind-strength compensate each other. We will come 
back to this finding and probable consequences in Sect. 17.51 
• The largest r e (j spread seen in our analysis occurs for the 

02 stars, with ranging from 44 to 55 kK when account- 
ing for all luminosity classes. Again, the more enriched of 
the two dwarfs (BI253 vs. BI237) is the hotter one, since 
the [N] effect dominates over the larger wind-strength. For 
the two supergiants (Sk-67° 22 and LH 10LW3-19), we 
find a similar T e ff difference, consistent with our predic- 
tions for a combination of [N] and log Q. Interestingly, the 
effect from a lower logg (3.70 vs. 3.90), as seen for the 

03 dwarfs, becomes inhibited by the extreme mass-loss. In 
Fig. [7] we noted that the N iv/N m emission line ratio be- 
gins to loose its sensitivity on \ogg at the 'F' -series of our 
model-grid. The wind-densities of the two supergiants are 
even higher (roughly corresponding to 'G'), and in this sit- 
uation the surface gravity does not play any role for deter- 
mining the line ratio. Finally, the difference between two 02 
giants (Nl 1-026 from the LMC and NGC 346-355 from the 
SMC), Ar e ff = 3 kK with respect to the 'average' solution 
for NGC 346-355, is larger than expected, since all other pa- 
rameters except for Z are quite similar, and the Z effect alone 
should amount to 1 kK (see Fig. [6j>. Note, however, that the 
errors in r e jf for both objects are larger than typical ( > 3 kK), 
related to the problems we encountered for ON2 (III) stars. 

The ON2 (III) stars. Let us point out already here that there 
seems to be a general, severe problem with the ON2 III class, 
as indicated from our results for Nl 1-031 and NGC 346-355, 
and the inspection of LH 64-16 (see Sect. 16.11 ). This problem 
must relate to some unknown physical process allowing for the 
presence of strong N in, N iv and N v lines in combination with 
weak, but still visible He i. Having explored hundreds of mod- 
els (both from our grid and additional, fine-tuned ones), vary- 
ing the N abundances as well as other parameters, it turned out 
that these lines could not be synthesized in parallel. Moreover, 
we explored a variety of potential sources also from the cmfgen 
side, manipulating the photosphere-wind transition zone, includ- 
ing wind-clumping with a variety of clumping laws, and near- 
photospheric X-ray s0 etc., an d found that none of these would 
solve the problem. Thus, we can be also sure that the problem 
is neither related to fastwind nor to our present nitrogen model 
atom. 

Also the mass-discrepancy found for these objects underpins 
their problematic nature. Using the hotter solutions for Nll- 
031 and NGC 346-355, we derived spectroscopic masses of 60 
and 53 Af , respectively. Ev olutionary masses from non-rotating 
tracks by ICharbonnel et al l (1 19931) and ISchaerer et ail (I1993I FI 
yield 90-95 M e for both stars. Note however that also for a cooler 
solution for NGC 346-355, at r eff = 49.5 kK, already Mas05 
stated quite a mass discrepancy, in this case 50 (spectroscopic) 
vs. 75 M Q (evolutionary). For another LMC ON2 III star, LH 64- 
16, not analyzed during the present work, Mas05 found an even 
larger discrepancy, 26 vs. 76 M . The authors suggested that the 
latter star might be the result of binary evolution, since it shows 
highly processed material at the surface, and since this star ap- 
peared to be located to the left of the ZAMS. This seems also 
to be the case for Nl 1-031 (at least for the hot solution), but not 

21 which do not help because of the destruction of Nm and Hei, see 
Paper II. 

22 For the hottest stars, the inclusion of rotation has only a modest ef- 
fect, resulting in ~ 10 % lower masses, e.g., Mas05, as long as the initial 
rotation is not high enough to induce quasi-homogeneous evolution. 
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for NGC 346-355. A similar problem was also found for som e 
of the close binaries in the R136 cluster (Mas sev et al.l 12002). 
supporting the idea that these stars suffered from binary inter- 
actions. An alternative to binarity might be quasi-homogeneous 
evolution, but the low v sin i values measured for our objects (on 
the order of 100 km s _1 ) render this possibility as unlikely. 

7.3. Niv/Nm emission line ratio 



To further test the significance of the Walbo rn et al.l |2002) 
classification scheme, we investigated the relation between the 
N iv/N in emission line ratio and T e s for our 02-04 sample stars 
in a quantitative way, similar to previous work by Mas05, but 
using our updated parameters (see Sect. 16.41 . Unlike Mas05 
who used EW ratios, we employed the line-strength ratios to 
be consistent with the classification scheme. Later on, we also 
consider an alternative line ratio, N V/14603/Niv/16380 (hereafter 
N v/N rv), to check its potential for classification purposes, and 
to test whether this ratio might break potential degeneracies in- 
herent to N iv/N m. 

Table [6] lists the line ratios for our targets, as derived both 
from the observations and from the synthetic models associated 
to our best-fitting solutions. Again, we discarded Nl 1-031 and 
R136-040 (see Sect. 17. 2b . Note that at the lower and upper end of 
the scheme the errors on the observed line ratios can become sig- 
nificant, due to noise and/or absence of N iv/14058 or N ni/14640. 
In those cases, we provide corresponding lower or upper limits 
and their uncertainties. Typical errors are on the order of 0.1 to 
0.2 dex. 

Figure [13] displays the derived effective temperatures as a 
function of the observed (upper panel) and 'model' (lower panel) 
N iv/N in line ratios, expressed logarithmically. Number designa- 
tions for each object refer to Table |6] Different colors are used 
for each spectral type: 02 (black), 03 (red), 03.5 (blue), and 
04 (purple). Symbols referring to luminosity classes and filling 
referring to [N] are as in Fig. [12] 

If we examine the observed emission line ratio, we find quite 
a monotonic behavior, confirming its T e s sensitivity, and similar 
trends as in the previous section. On the one hand, dwarfs and 
giants behave rather similar (though the giants seem to be a lit- 
tle cooler), except for the 02 types where the spread is larger. 
On the other, the observed relation for supergiants lies in paral- 
lel to the relation for dwarfs, but at considerably lower r e ff . As 
already pointed out, this is caused by the different gravities and 
wind-strengths. E.g., from Table [5] we find log Q values in the 
ranges [-13.00, -12.53] for the LMC dwarfs, around -13.2 for 
the SMC dwarfs, [-12.76, -12.36] for the giants, partly overlap- 
ping with the dwarfs, and [-1 1.97, -1 1.60] for the supergiants. 

The various spectral subtypes are locate d in quite different 
ranges (by definition, consistent with the Wal born et al.l |200"2| 
scheme): All 02 stars are located at log N iv/N in > 0.4, and 03 
dwarfs and 03.5 giants around log Niv/Niii ~ (Niv ~ Nm). 
Finally, one of the 04 dwarfs is located at log N iv/N in « -0.7, 
for the other one (#8) we are only able to provide an upper limit, 
log N iv/Nm £ -0.3, whilst the 04 supergiant (#15) lies close to 
the 03.5 border. 

We see small [N] effectS on the line ratio, where a larger 
abundance tends (but not necessarily) to increase the emission 
line ratio within a given spectral subtype, in particular for the fol- 
lowing pairs of stars: 03.5 dwarfs - LH 101:W3-24 (#6, upper 
limit) and LH 8LW28-23 (#5); 03 dwarfs - AV435 (#4, SMC) 
and Nl 1-060 (#3); 02 dwarfs - BI237 (#2) and BI253 (#1); 02 




log (NIV4058/NIII4640)„ 



Fig. 13. Effective temperatures as a function of log N iv/N m, for 
02-04 sample stars. Spectral types distinguished according to 
legend. Upper panel: observed emission line ratios, limits indi- 
cated by arrows; lower panel: emission line ratios predicted from 
the (globally) best-fitting models. Objects numbered according 
to Table [6] Filling and symbols as in Fig. [12] Dashed vertical 
lines (at log Niv/Niii = -1, 0, 1) correspond to the limits used 
in Fig. [6] For object #9 (NGC 346-355), r eff and theoretical line 
ratio according to our 'average' solution. Stars #4, #8, and #9 are 
SMC stars. 



supergiants - LH 10LW3-19 (#14) and Sk-67° 22 (#13). Since, 
on the other hand, the emission line ratio should decrease for en- 
hanced [N] when keeping all other parameters fixed, a consistent 
solution requires these objects to have a higher T e ff. For the pair 
of 04 dwarfs - AV177 (#8, SMC) and LH 8LW28-5 (#7), we 
note that the star with the higher abundance (#7) seems to have a 
lower emission line ratio.Pl but this is consistent with the similar 
T e ff of these stars. 

Indeed, there are not only [N] effects, but also mass-loss ef- 
fects on the emission line ratio. For all pairs considered above, 
the object with higher [N] has also a higher wind strength, 
log Q (cf. Table [5), which potentially counteracts the [N] effect. 
Comparing the differences in [N] with the differences in log Q, it 
turns out that the [N] effect should dominate in all cases though. 

Taken together, our findings explain the almost monotonic 
increase of T e g with N iv/N m also within the individual spectral 
subtypes, consistent with the derived scatter of r e ff per subtype. 



slightly contaminated by mass-loss effects, see below. 



but note also the quite large uncertainty for object #8. 
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Table 6. Observed and predicted line-strength ratios for N iv/N m and N v/N iv, for the 02-04 stars analyzed within this work. 
Observed ratios (or limits) inclusive errors. Numbers refer to Figs. [T3landfl4l 



Star 


# 


SpT 




log( 


NIVn 
NII1- 1 


i / N 1 V \ 




log( 


NV x 
NIV ' 


i / NV \ 


Comments 








(kK) 


obs 


model 




obs 


model 




BI253 


1 


02 V((f* )) 


54.8 


0.67 


±0.15 


1.85 




0.53 


±0.23 


0.43 


N ni/14640 and NiV/14058 underpredicted 


BI237 


2 


02 V((f )) 


53.2 


0.59 


±0.14 


1.58 




0.36 


±0.14 


0.21 


N ni/14640 underpredicted" 


Nl 1-060 


3 


03 V((f )) 


48.0 


0.08 


±0.02 


-0.03 




-0.07 


±0.03 


-0.30 


Nrv/14058 and N V/14603 underpredicted' 1 ''' 


AV435 


4 


03 V((f )) 


46.0 


0.01 


±0.18 


-0.03 


< 


-0.30 


(+0.10) 


-0.46 


N V/14603 diluted in noise 


LH 81:W28-23 


5 


03.5 V((f+)) 


47.0 


-0.16 


±0.04 


0.09 




-0.01 


±0.05 


-0.18 


N ni/14640 underpredicted 


T T T 1 Al \T T~) >~\ A 

LH 101:W3-24 


6 


03.5 V((f+)) 


47.0 


< -0.20 


(+0.12) 


-0.12 


< 


-0.28 


(+0.14) 


-0.38 


N in/14640 underpredicted 


LH81:W28-5 


7 


04V((f+)) 


44.0 


-0.66 


±0.10 


-0.55 




-0.53 


±0.13 


-0.60 


Satisfactory fits 


AV 177 


8 


04 V((f)) 


44.0 


< -0.26 


(+0.13) 


-0.34 


< 


-0.36 


(+0.09) 


-0.57 


Only N ni/14640 and N iV/16380 visible 


NGC 346-355 


9 


ON2 III(f ) 


53.0 


0.91 


±0.09 


1.33 




0.19 


±0.04 


0.21 


NiH/14640, NiV/14058, Nv/14603 underpred/ 


Nl 1-026 


10 


02 III(f ) 


49.0 


0.54 


±0.06 


0.73 




0.14 


±0.02 


-0.08 


N ni/14640, N iv/14058, N v/14603 underpred."" 


R136-018 


11 


03 ni(F) 


47.0 


> +0.22 


(-0.21) 


-0.18 


> 


-0.20 


(-0.13) 


-0.32 


Nni/14640 and NiV/14058 underpredicted 


LH 90: ST 2-22 


12 


03.5 III(f+) 


44.0 


-0.18 


±0.06 


-0.13 




-0.20 


±0.04 


-0.39 


Satisfactory fits 


Sk-67° 22 


13 


02 If /WN5 


46.0 


0.87 


±0.14 


1.04 




0.35 


±0.03 


0.31 


N V/14603 under, NiV/14058 overpredicted 


LH 10LW3-19 


14 


02 If 


44.0 


0.40 


±0.07 


0.41 




0.05 


±0.14 


0.00 


N v/14603 overpredicted 


Sk-65° 47 


15 


04 If 


40.5 


-0.19 


±0.12 


0.23 


< 


-0.38 


(+0.08) 


-0.52 


N iV/14058 overpredicted 



(fl) For fits, see Paper II, Appendix C; Compromise solution, see Paper II; <c,) 'Average' solution. 

Notes. Errors in brackets provide uncertainties of lower or upper limits. Predicted line ratios ('model') drawn from best-fitting synthetic spectra. 
Nl 1-031 and R 136-040 discarded, see text. 



Differences in background metallicity and wind-strength seem to 
play a secondary role, compared with the larger impact of [N] . 

If we now inspect the line ratios predicted by our best-fitting 
models (lower panel), we mostly find quite similar values and 
trends. For the bulk of the stars there are only small shifts due 
to minor problems in properly fitting the lines (see Table |6] for 
particular comments on each star). However, more severe differ- 
ences are present for BI237 (#1), BI253 (#2), and NGC 346-355 
(#9, SMC). For the former 02 dwarfs, we predict too weak N in 
emission, thus overestimating the line ratio, but note that the ob- 
served emission is also rather weak for these stars. For NGC 
346-355, the problem is different since we are not able to re- 
produce both lines using the 'average' solution (see Fig ]B.12l > 
displayed here. 

We conclude that, for a given luminosity class, the effective 
temperatures are a rather monotonic function of log N iv/N m, 
that the scatter within a spectral subtype is mostly due to abun- 
dance effects, and that our models are in fair agreement with the 
observed line-ratios, except for the hottest objects where we un- 
derestimate the observed (low) N m emission strength. 

7.4. N v/N iv line ratio 

Since both NiV/16380 and Nv/i4603 are absorption lines (i.e, 
less affected by complex formation processes), since they turned 
out to be quite reliable during our analyses, and since N v/14603 
is very r e jf sensitive (Sect. [5), we checked the corresponding line 
ratio as a potential diagnostic tool, which might be even used for 
future classification purposes. 

From Fig. [14] we see that the relation T e ff vs. log Nv/Niv 
is remarkably monotonic. By inspection of the observed line 
ratios (upper panel), we find again two different trends, one 
for dwarfs and (bright) giants and another one for supergiants. 
The objects are basically grouped together within three regions: 
02 stars with log (Nv/Niv) obs > 0, 03/03.5 dwarfs and gi- 
ants around log (N v/Niv) b s ~ -0.3. . .0, and the 04 stars with 
log (N v/N iv) bs -0.4. The only discrepant object seems to 
be the SMC 03-dwarf AV 435 (#4), which appears at the edge 



60 



50 



45 



40 



. □ 1 

. Alll 
. OV 


1 1 1 1 1 










1 A9 


01 

02 






! A10 

45 


■ 13 




•7 — 8 A12 








■ D15 

1 




02 - 

03 " 

03.5 - 

04 



-0.5 0.0 0.5 

log (NV4603/NIV6380) 01 „ 



60 



50 



40 



. □ I 


1 1 I 1 1 1 1 I 1 




. A III 






. OV 










Ol 






9A>2 




A10| 






•3 






06A11 «5 . 






; 


■ 13 




7 8 A12 UH 






□ 15 | 






02 - 






03 " 






03.5 " 




1 .... 1 


04 




-0.5 0.0 


0.5 




log (NV4603/NIV6380) mo 


dels 



Fig. 14. As Fig. Q/5] but for the N v,!4603/NiV/t6380 line ratio. 



of the 04-region. However, this 'erroneous' position could be 
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tracked down to a considerable error in the measured line ra- 
tio, because of very weak Nv lines diluted in the continuum 

(Fig. EES. 

Compared to the N iv/N m emission line ratio, there seems 
to be a clearer separation between the different subtypes, e.g., 
Sk-65° 47 (#15) is located closer to the remaining 04 ob- 
jects, whilst regarding N iv/Nm it is closer to the 03 V/03.5 III 
group. Moreover, there is a clear separation between the two 02 
dwarfs (#1,2) and the ON2 giant (#9, SMC) because of weaker 
Niv/16380. 

Comparing now with the predicted line ratios (lower panel), 
we see that predicted and observed ones agree quite well over 
the complete range. Still, there are certain shifts because of a 
non-perfect representation by our synthetic lines (Table |6j, but 
interestingly we do no longer find the extreme differences for 
BI237, BI253, and NGC 346-355 as present in the Niv/Nm di- 
agram, which indicates that we are able to obtain a satisfactory 
representation of the N iv/N v ionization balance. So far, we can- 
not comment on any [N] related bias, since we did not perform 
corresponding theoretical studies. In case the N v/N iv line ratio 
might be used in future classification schemes, this will become 
certainly necessary. 

The monotonic behavior of the N v/N iv line ratio illuminates 
its promising potential, particularly for the hottest objects (02), 
simply because there is more N v than N m present at these tem- 
peratures. We suggest to investigate this possibility when larger 
samples become available. 

7.5. Caveats for low-metallicity stars 

Summarizing our fin dings, we conclude that already the present 
Walb orn et al.l (120021) classification scheme allows for a reason- 
able relation between spectral type and effective temperature, as 
long as it possible to discriminate the luminosity class. The only 
significant bias in the scheme might be produced by nitrogen 
abundance effects (or extreme variations in wind strength). E.g., 
if the nitrogen abundance is not the same in an 03 III and an 
03.5 III star, then the 03 III star is not necessarily hotter than 
the 03.5 IE. 

However, there are also important caveats regarding low- 
metallicity (e.g., SMC) stars. As it is well-known, effective 
temperatures increase with decreasing Z for a given spectral 
type if the classification is based on the helium ionization bal- 
ance, i.e., the Hei/ He n line-st rength ratio (04-09.7 stars, e.g., 
Bouret et aj] |2003l Mas04/05. [Heap et alj 120061 iMokiemetai] 
2006, 2007b). According to our predictions, this no longer needs 
to be true for the earliest O-stars classified by means of nitrogen. 

(i) 02/03 stars. Though for similar [N] a lower metal- 
licity implies a higher T e g (at least for spectral types 
03(dwarfs)/03.5(giants and supergiants) and earlier, see 
Fig- 01, this effect might be counteracted by a different nitro- 
gen content if we assume typical (maximum) enrichments, 
of roughly +0.6, +0.9, and +1.0 dex above the correspond- 
ing MW^J-iMC, anclSMC baseline abundance, as predicted 
by iBrott et all (1201 lab for a 40 M Q star at an initial rota- 
tion of 270 kms~'. As visible from Fig. [8] for an 03 dwarf 
with log g — 4.0 dex, our predictions indicate T e ff » 47-48 kK 
for the Galaxy and r eff * 46 kK for the LMC/SMC0 Such 
cooler or at least similar r e jj for stars in a lower Z environ- 
ment should be present only for typical nitrogen abundances; 



25 See also Fig. [12] where the SMC 03 dwarf AV 435 indeed is cooler 
than the corresponding LMC objects. 



objects with a considerably different enrichment will con- 
tribute to enlarging the spread, 
(ii) 03.5/04 stars. In view of the results from Sect. 14.2.11 
there might be an additional problem around 03.5/04. From 
Fig. [6] SMC stars should be cooler than corresponding LMC 
ones, even at a similar [N], for a line ratio around N iv/N m 
0. 1 (which represents a lower limit for 03.5 V and 04 I/III). 
Since the SpT-r e fif relation is rather monotonic for LMC stars 
at all spectral types, derived either by helium or nitrogen 
(see Fig.[T2l. the corresponding relation for SMC stars might 
be not monotonically connected between these two regime: 
Even though there seems to be a monotonic relation between 
N iv/N in and T e ff on the hot side, the potential 'jump' would 
be caused by the fact that the SMC SpT-r e ff relation on the 
cooler side (based on He i/He n) lies above the LMC relation. 

To confirm or disprove these predictions and caveats, a thorough 
analysis of a large sample of SMC stars is certainly required, 
given the few SMC objects investigated so far. 

7.6. SMC Of -stars 

The luminosity criteria for O-type stars are based primarily 
on their "Of" characteristics; i.e., the strength of the Nm 
triplet and the He ii/14686 emission^ As already emphasized by 
Mas04/05/09, this can lead to significant problems with lumi- 
nosity classification of O-stars found in low metallicity environ- 
ments, such as the SMC. 

The problem is that the emission strengths of the Of features 
are tied to mass-loss rates (both for Hen and Nm), which in 
turn scale with metallicity, and to (average) nitrogen abundances 
(for N m)0One thus expects a star in the SMC to show weaker 
Of-characteristics for the same physical properties than would a 
similar star in the Milky Way. As demonstrated in Sect. 13.11 in 
particular the N in emission strength of the SMC star should be 
lower, unless the nitrogen abundances of the two objects were 
similar, which would mean the rare case of an extreme enrich- 
ment of the SMC star. 

Throughout their studies of Magellanic Cloud O-type stars, 
and particularly those in the SMC, Mas04/05/09 found numer- 
ous examples where the Of-type properties indeed were weaker 
than would be expected given the absolute magnitude of those 
starsF^I. We can demonstrate the weakness of Of-features sta- 
tistically as follows. If we restrict ourselves to the SMC O- 
type stars with spect ral subtypes d etermined from slit spec- 
troscopy (Table 6 of iMassevf 120021) we find 74 stars, only 5 
(7%) have "f " type designation s . Sim ilarly, of the 83 O-type 
stars listed by lEvans & Howarthl (120081) . only 7 (8%) have any 
"f" designation, and all of these are either "(f)" or "((f))" in- 
dicating that He H/14686 is weakly in emission or in absorp- 
tion. In contrast, Of-characteristics abound among O-type stars 
in the M ilky Way. Of the 378 Gal actic O-type stars cata- 
logued bv lMafz-Apellaniz et al.l ( 120041) . 160 (42%) display Of- 
characteristics. Since Of-type stars are brighter than non-Of stars 
(at least in the Milky Way), this may overestimate the true per- 

26 These emission fea tures were link ed with luminosity for Galactic 
stars (see, for example, Walborn 1972), a refinement over previous lu- 
minosity criteria based primarily on the SiiV/14089 to He 1/14143 ratio. 

27 Counteracting the increase of N in emission because of less block- 
ing in a low Z environment. 

28 The alternative explanation that all such discrepant stars were too 
bright owing to their being binaries was contradicted by the excellent 
fits obtained to the spectral features; this would require both compo- 
nents of a binary to be of identical spectral subtype and brightness. 
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centage. If we instead restrict ourselves just to the sample of 
24 O-type stars with well- established distances by Hipparcos 
(Maf z-Apellaniz et al.ll2004l Table 7) we find that 7 (29%) have 
Of-characteristics. 

Thus we conclude that there are strong indications for a 
significantly lower percentage of SMC Of-stars, compared to 
Galactic ones, which we attribute to the effects outlined in 
Sect. 13.11 Remember, however, that a lower wind-strength alone 
is not sufficient to explain this finding (Fig. [2] upper panel), but 
that a lower average nitrogen content needs to be present as well. 

8. Summary and conclusions 

We investigated open questions raised by our previous studies 
on the formation of N ni/14640 and N iV/14058. We provided first 
theoretical predictions for the N iv/N m emission line ratio, and 
confronted these predictions with observational findings, con- 
centrating on a sample of early-type O-stars. The results of this 
work can be summarized as follows. 

1. The emission strength of the Nm triplet from objects with 
similar T e ff and log g depends on their metallicity, associated 
mass-loss, and nitrogen content. Whilst even under SMC 
conditions lower mass-loss rates alone are not able to com- 
pensate the increase in emission for decreasing Z, a lower 
[N], coupled to a significantly lower base-line abundance, 
can easily outweigh the Z effect and lead to overall lower 
emission strengths. This might explain the relatively low 
number of SMC 'Of stars. 

2. Our models predict an only weak Z-dependence of 
N iV/14058 (contrasted to the N in triplet). SMC-abundance 
models with T e g < 45 kK display slightly more emission than 
their Galactic counterparts, and vice versa for hotter temper- 
atures. Much stronger is the impact of wind-strength though. 

3. It turned out that Niv/14058 behaves quite unexpectedly 
when [N] is increased in low-M models. For almost the 
whole temperature range, we either obtain more absorption 
or less emission, compared to models with a lower nitro- 
gen content. For a specific temperature range, N iV/14058 can 
even switch from emission to absorption when increasing 
[N]. For Galactic stars at 44 kK < r eff < 50 kK and com- 
paratively low M, our models imply that if N iV/14058 is ob- 
served in absorption, this would indicate a strong nitrogen 
enrichment. 

4. We provided first theoretical predictions on the N iv/N m 
emission line ratio, as a function of Z, log Q, and [N], by 
studying line-ratio iso-contours in the r e ff-log g plane. For an 
emission line ratio of unity (i.e., a spectral type of 03.5 I/III 
or 03 V), the corresponding increases with [N] (~1 kK 
per increment of 0.2 dex in [N]) and log g (~1 kK per incre- 
ment of 0. 1 dex in log g). In addition, it should decrease with 
Z, at least for higher T e g (~2-3 kK difference between SMC 
and MW objects), and log Q (~2-4.5 kK between low- and 
high-M models). 

5. We performed a comparison with results from the alternative 
model atmosphere code cmfgen, for a small grid of early O- 
type dwarfs and supergiants. Our basic predictions regard- 
ing the impact of [N] on (i) N iv/i4058 for low-M models 
(see item 3), and (ii) on the N iv/N in emission line ratio (see 
item 4) were confirmed by corresponding cmfgen results. 
Regarding specific line predictions, we found a mostly sat- 
isfactory agreement, except for some systematic deviations: 
For early O-stars, fastwind produces more emission at the 
Nm triplet, less emission at Niv/14058, and mostly much 



more absorption at N iV/16380. This would lead to lower T e g 
and quite different [N] in analyses performed by means of 
cmfgen, if concentrating on the Nm triplet and Niv/14058 
alone. Fortunately, the remarkably good agreement of the 
H/He lines in both codes enables an identification of poten- 
tial problems regarding the nitrogen lines, as long as the he- 
lium ionization balance is used to constrain T e ff . In the hot O- 
star domain, the latter approach is no longer feasible because 
of vanishing Hei. Nevertheless, potential problems should 
become obvious also here when relying on the N V/M4603- 
4619 doublet, which turned out to be very sensitive on T^g 
as well as code-independent. 

6. We confronted our theoretical predictions with results from 
an analysis of a medium-size sample of LMC/SMC O-stars, 
drawn from studies by Massey et al. (early types), and from 
Paper II. The basic difference to the Massey et al. analyses 
is found in the procedure for deriving T e ff, where we used 
the nitrogen ionization balance in the hotter T e g regime in- 
stead of the helium one, to avoid any degeneracy. For the 
cooler objects of our sample (r e ff < 44 kK) we mainly relied, 
when possible, on helium, using nitrogen as a consistency 
check. For these stars we found similar or slightly cooler T e g 
compared to Massey et al.. Considerably hotter r e ff, on the 
other hand, were inferred for the earliest O-stars, by means 
of the nitrogen diagnostics. Nevertheless, in most cases the 
corresponding synthetic He n/He i lines were still consistent 
with the observations, or indicated only slightly lower tem- 
peratures. Notable exceptions are the ON2 III stars (see be- 
low). The nitrogen abundances derived within our analysis 
are consistent with our results from Paper II: again, the bulk 
of the stars displays a considerable enrichment. 

7. By inspecting the inferred effective temperatures, we saw 
that r e ff increases from supergiants to dwarfs for all spectral 
types, consistent with earlier results. For spectral types later 
than 03.5 (down to 09.5), LMC giants are cooler by ~1 kK, 
and supergiants are cooler by ~4 kK, compared to dwarfs. 
For types earlier than 03.5, this difference (and also the scat- 
ter) becomes larger, amounting to ~10 kK at 02 when com- 
paring supergiants and dwarfs. For LMC dwarfs and giants 
later than 03.5, and for all LMC supergiants, we found lin- 
ear relations between and spectral type, again consis- 
tent with previous work. The earliest dwarfs and giants, on 
the other hand, display a much steeper increase in r e jf . The 
dominating effect responsible for the scatter in the SpT-r e ff 
relation at earliest types was attributed to difference in [N], 
where for a given spectral type more enriched objects are 
typically hotter. 

8. The relation between the observed Niv/Nm emission line 
ratio and r e jf turned out to be quite monotonic, if discrim- 
inating for luminosity class. Because of the high T e g de- 
rived for the earliest stars by means of the nitrogen ioniza- 
tion balance, we did not find the pronounced degeneracy of 
the Niv/Niii emission line ratio as claimed by Mas05. The 
scatter found within a spectral subtype is, again, primarily 
produced by abundance effects. Our model predictions are in 
fair agreement with the observed line-ratios, except for the 
hottest objects where we underestimate the observed (low) 
Nm emission-strength. 

9. We provided first insights into the relation between T e ff and 
the N V/L14603-4619/N iv/16380 absorption line ratio, which 
is remarkably monotonic, particularly for the hottest objects 
in our sample, and we highlighted the promising potential of 
this line ratio for future classification schemes. 
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Both our theoret ical prediction s and o ur observational analysis 
suggest that the Walbor n et al ] (120021) classification scheme is 
able to provide a meaningful relation between spectral type and 
effective temperature. In particular, and as one might have ex- 
pected, the N iv/N m emission line ratio changes with effective 
temperature, all other factors being equal. However, this ratio 
is also sensitive to surface gravity, mass-loss rate, and to ni- 
trogen abundance, which are expected to vary among a sam- 
ple of stars. Thus, the significance of the classification scheme 
(within the uncertainties caused by nitrogen abundance) might 
be only warranted as long as it is possible to fairly discrimi- 
nate the luminosity class (as a proxy to gravity), and as long 
as there is a strong correlation between spectral type/luminosity 
class and wind-strength. If, e.g., there would be weak-winded 
stars dMarcolino et al.ll2.009L iNaiarro et al.ll2.Ql II and references 
therein) also in the early O-type regime and not only at later 
spectral types, the monotonicity with respect to T e ff might be- 
come severely disturbed. 

A clear identification of early O-type luminosity classes from 
spectral morphology alone becomes difficult in low-Z environ- 
ments such as the SMC. We emphasize the same point as made 
by Mas04/05/09, that the standard luminosity classification crite- 
rion, primarily based on the morphology of He H/14686, is signif- 
icantly biased on mass-loss rates. Owing to lower wind-strengths 
for SMC conditions, Heii/14686 is typically in absorption not 
only for lc V, but also for lower gravity objects which under 
Galactic conditions would correspond to lc I/III. Thus, there are 
O-type stars in the SMC whose physical properties (visual lu- 
minosities and surface gravities) might be in accord with them 
being giants or supergiants whereas their spectroscopically de- 
termined luminosity classes may be dwarfs or giants, respec- 
tively. To circumvent this caveat, other information, such as the 
absolute visual magnitude of the system, may have to be ap- 
pealed to. Without such additional information, an appropriate 
classification may rest on a more precise determination of the 
surface gravity, e.g., using a visual inspection of the wings of 
H r . Another possibility is to extend the classification scheme by 
including the strength of H„ (which needs to be carefully cali- 
brated, since also this wind-line remains in absorption). 

Our study also implies important consequences which need to 
be validated or investigated in future work: 

(i) The SpT-r e ff scale constrained from our LMC sample turned 
out to be, for a given luminosity class, more or less mono- 
tonic. However, the majority of the analyzed objects dis- 
played a considerable nitrogen enrichment. To quantify the 
actual impact of nitrogen abundance would require the anal- 
ysis of a significant number of un- or mildly enriched earliest 
MW/LMC objects^ that according to our predictions should 
be cooler than the enriched ones. 

(ii) Typical SMC stars of earliest spectral types might have ef- 
fective temperatures below corresponding LMC objects, and 
the overall SpT-r e fif relation for SMC stars might be non- 
monotonic around 03.5/04. 

(iii) Our predictions suggest that there should be (enriched) 
Galactic 04 dwarfs with emission at Nni/14640 and ab- 
sorption at Niv/140 58 (item 3 from above). So far, the 
IWalborn et alj (120021) classification scheme only states that 
Niv/14058 should be absent for such objects. A quick 

29 For SMC objects, the nitrogen baseline abundance is too low to al- 
low for a clear-cut classification and spectroscopic analysis by nitrogen 
lines, at least in the optical. This was already a problem for unenriched 
LMC stars, e.g., R136-040, and the similar objects R136-033 and R136- 
055. 



inspection of the IACO B database of Galactic OB stars 
( Simo n-Diaz et al.ll201 ll) allowed us to find a first indication 
of this effect, for HD 46223 (04 V((f))), and we were able 
to reproduce its strong emission at N ni/14640 together with 
a pronounced absorption at N rwi4058, for a considerable ni- 
trogen enrichment, 
(iv) As became apparent throughout this work (first indications 
were already found in Paper II), we encountered severe prob- 
lems for ON2 III stars. For these objects, we were not able 
to find a simultaneous fit of the pronounced N m/N iv/N v 
lines and the weak, but still clearly present Hei/l 4471, in- 
dependent whether we used fastwind or cmfgen. Though we 
were able to derive a cooler (fitting N m/N iv and He i) and 
a hotter solution (fitting N iv/N v), we were not able to fa- 
vor one of those. In some cases, the T e ff difference between 
the hotter and the cooler solution is extreme, e.g., for Nll- 
03 1 analyzed in Paper II. Because of the restricted quality of 
our present dataset, new observations with higher S/N and 
higher resolution would be extremely valuable to verify our 
results using the nitrogen ionization balance and/or the pure 
He n/He i results from Massey et al. when He i (and N m) be- 
come extremely weak. 

Future work is certainly needed to address the aforementioned 
issues. Upcoming analyses of extensive O-star samples, e.g., 
from the VLT-FLAMES Tarantula survey and the IACOB data 
base, will be fundamental for trying to explain these open ques- 
tions. 
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Appendix A: Comparison with cmfgen 

Figures lATI to IAT81 provide a detailed comparison between H/He and Nm/Niv/Nv cmfgen spectra for models d2v, d4v, s2a, 
and s4a (see TableO, and corresponding fastwind profiles from closest or almost closest grid models, for a nitrogen abundance of 
[N] = 8.78 and [N] = 7.78. If not explicitly stated else, no convolution has been applied to the spectra. For details, see Sect. [5] 
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Fig.A.l. Model d4v at [N] =8.78. Comparison of H/He/N spectra from cmfgen (green) and fastwind, at the closest grid-model 
(black: T eE = 41 kK, logg = 4.0, log Q = -12.8, [N] = 8.78) and at neighboring grid models with r eff = 40 kK (red) and T eS = 
42 kK (blue). To allow for an easier comparison, all profiles have been convolved with v sin i = 30 km s~' . 



J.G. Rivero Gonzalez et al.: The earliest O-stars, Online Material p 2 



NIII4003 



NIII4097 



NNI4195 




4001.5 4002.8 4004.2 
wavelength (A) 



NIII4634 -4640-4642 




4630.0 4634.7 4639.3 
wavelength (A) 



1.03 
1.02 
1.01 
1.00 
1 0.99 
0.9B 
0.97 
4053.5 



NIV405B 





\ 


\ 


/ 



4056.2 4058.8 
wavelength (A) 



1.1 

£ 1-0 

Z 0.9 
a. 

c 0.8 

01 
OA 

i 0.7 
E 

* 0.6 

0.5 









\ / 



4095.3 4096.6 409B.0 
wavelength (A) 



NIII4510-14-18 



4507.0 4511.7 4516.3 
wavelength (A) 

NIV6380 




4193.7 4195.0 4196.4 
wavelength (A) 

NIV3480 








* 


n 

: 



3474.0 3478.7 3483.3 
wavelength (A) 

NV4603-4619 



6377.0 6379.7 6382.3 
wavelength (A) 



4601.7 4604.3 
wavelength (A) 




4607.0 



Fig. A.2. Model d4v at [N] =7.78 (solar). Comparison of N spectra from cmfgen (green) and fastwind, at the closest grid-model, 
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Fig. A.3. Model d2v at [N] =7.78 (solar). Comparison of N spectra from cmfgen (green) and fastwind, at the closest grid-model 
(black: r eff = 46 kK, logg = 4.0, log Q = -12.45, [N] = 7.78), and at the neighboring grid model with r eff = 47 kK (red). The H/He 
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Fig.A.4. Model d2v at [N] =8.78. Comparison of H/He/N spectra from cmfgen (green) and fastwind, at the closest grid-model 
(black: r eff = 46 kK, logg = 4.0, log Q = -12.45, [N] = 8.78) and at neighboring grid model with r eff = 47 kK (red). No convolution 
has been applied. 
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Fig. A.5. Model s4a at [N] =8.78. Comparison of H/He/N spectra from cmfgen (green) and fastwind, at the closest grid-model 
(black: r eff = 39 kK, logg = 3.5, log Q = -12.10, [N] = 8.78). 
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Fig. A.6. Model s4a at [N] =7.78 (solar). Comparison of N spectra from cmfgen (green) and fastwind, at the closest grid-model 
(black: T eS = 39 kK, logg = 3.5, log Q = -12.10, [N] = 7.78). The H/He spectra remain as in Fig. [A3] 
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Appendix B: Comments on the individual objects 

In the following, we give specific comments on the individual 
objects, regarding peculiarities and problems found during our 
analysis. We separate between galaxy membership, and sort by 
luminosity class and spectral type. Line fits are displayed in 
Figs. IO to EE! for important H/He/N lines: H a , % H r , H a , 
H e , UeiAA 4387, 4471, 4713, Hen(+Hei)^4026, HenXi4200, 
4541, 4686, 6406, 6527, 6683, NmAA 4003, 4097, 4195, 4379, 
AA4634 - 4640 - 4642, and AA45 10 - 45 14 - 45 1 8, N ivAA 4058, 
6380, and N V/L14603/4619. All spectra were corrected for ra- 
dial velocity shifts. 

If not explicitly stated, any comparison made in the follow- 
ing text refers to the results from Mas05. 

B.1. LMC stars 

R136-040 - 02-3.5 V (Fig.Ej}. Th is star could not been clas- 
sified by Mas05 using the scheme by Walb orn et al.l (2002), be- 
cause neither Nni/M4634 - 4640 - 4642 nor NiV/14058 were 
visible in the spectra. As outlined in Sect. 16.21 for the R136 stars 
we have spectra from both STIS and FOS available. Taking ad- 
vantage of the better quality of the STIS data and using H a , H r , 
and He n/14541, we derived a similar lower limit (no He i and no 
nitrogen lines)! on T e s as Mas04, but a substantially larger (by 
0.2 dex) gravity, which agrees better with its dwarf designation. 
Our analysis also resulted in a low helium content, Yu e = 0.08 
(Mas04: Yu e - 0. 1). Note the discrepancy in the cores of H<s, H e , 
and He H/l/14026, 4200 when comparing with the FOS data. Such 
a discrepancy was also found for the remaining H/He n lines in 
the FOS data, when used instead of the STIS spectra. 

The missing nitrogen lines imply an upper limit for the ni- 
trogen content corresponding to the LMC baseline abundance 
[N] = 6.90, which also agrees quite well with the low He con- 
tent. 

LH 81:W28-23 - 03.5 V((f+)) (Fig. |53}. The modeling of 
this star was straightforward, and we obtained similar results 
as Mas05. However, we considered a larger /? = 1.0 to better 
reproduce the marginal P-Cygni profile at He n/14686, which 
migh t indicate a luminosity class III object (see IWalborn et al] 
2002). To preserve the fit of H a , we needed to reduce M. All 
nitrogen lines are consistent with the temperature derived from 
the helium ionization equilibrium. The quite large nitrogen 
abundance ([N] = 8.40) agrees well with the helium abundance 
Y\ie = 0.25, indicating an evolved nature of this object. 



LH 10LW3-24 - 03.5 V((f+)) (Fig.HO). We derived a some- 
what cooler T e ff (by 1 kK) together with a lower helium content, 
^He = 0.10, which was consistent for all helium lines. Nitrogen 
lines are barely visible, because the spectrum of this star displays 
more noise than the bulk of our sample stars, caused by the use 
of a narrow extraction aperture to reduce effects from nebular 
emission for the ground-based observations (cf. Sect 16.2) . Due 
to the noisy spectrum, we were only able to infer an upper limit 
for the nitrogen abundance, [N] < 7.78. 

LH 81:W28-5 -Q4 V((f+)) (Fig.|E!>. This is one of the stan- 
dards used bv lWalbornetal.l (120021) for defining the 04 V((f+)) 
class. A consistent analysis of the helium and nitrogen ioniza- 
tion equilibrium yielded a cooler temperature, T e g = 44 kK, 



which required a helium abundance of Yn e = 0.15 to reproduce 
HeI/14471. An excellent fit to most nitrogen lines from all three 
ionization stages was achieved for this star, indicating a signifi- 
cant enrichment, [N] = 8.38. 

R136-018 - 03 III(f ) (Fig. |B3]>. Also for this 03 giant we 
have used data from both STIS and FOS. A consistent analysis 
of the H/He lines from STIS and nitrogen lines (mostly from 
FOS) suggested a hotter T e g, by 2 kK, as well as a higher surface 
gravity, by ~0. 1 dex. Again, we found discrepancies in the cores 
of the H/He n lines from the FOS spectra, except for He H/14686. 
An acceptable fit for N m/N iv/N v using [N] = 8.18 was possible, 
where only N iV/14058 was slightly underpredicted. Even though 
the nitrogen lines contained in the STIS dataset, N 111/1/145 10 - 
4514 - 4518 and Niv/16380, are diluted in the continuum, they 
support our analysis. 

LH 90:ST 2-22 - 03.5 III(f+) (Fig. HTSI An unproblematic 
analysis provided the same results as obtained by Mas05, except 
that we opted for a lower helium abundance, Yu e = 0. 15. Again, 
a remarkable fit to the nitrogen lines was possible, at [N] = 8.58, 
indicating an extreme enrichment. 

Sk-67° 22 - 02 IF/WN5 (Fig.lBJl This star was re- classi- 
fied3 as 02 If/WN 5 by ICrowther & Walbornl d20lT1) using 
their updated classification scheme, because of the Hp P-Cygni 
profile. Using lines from N m/N iv/N v, we inferred T e g = 46 kK 
which is hotter than the lower limit (from very weak He 1/14471) 
quoted by Mas05. Our fit seems to slightly overpredict the emis- 
sion in N iV/14058 and to underpredict the N v doublet. An ex- 
treme nitrogen abundance, [N] = 8.78, was required, the largest 
one found in our sample. Such an abundance would be certainly 
too large when comparing even with strongly nitrogen-enhanced 
O-stars, and also with predictions from e volutionary calculations 
tailored for the LMC (IBrott et al.ll201 1 bl and paper II), thus sup- 
porting a rather evolved nature of this object and its 'slash-star' 
design ation. This star was also analyzed by iDoran & Crowtherl 
d201 ll) using N rv/N v lines (without discussion of He i and Nm), 
only providing a T e g = 49.3 kK for this object. Such hotter tem- 
perature would improve our fits for N iV/14058 and the N v dou- 
blet, but is inconsistent with the observed strength of N m. 

LH 10LW3-19 - 02 If (Fig. EE). For this supergiant, a con- 
sistent He/N analysis allowed us to derive T e g = 44 kK, hotter 
than the lower limit (marginal HeI/14471) assigned by Mas05. 
Using N m/N iv/N v in parallel, we achieved an almost excellent 
fit for the nitrogen lines at [N] = 8.18. 

Sk-65° 47 - 04 If (Fig. HOI The parameter set derived for 
this star using H/He/N lines is quite similar to the results from 
Mas05, with somewhat larger Yn e = 0.12, A potential discrep- 
ancy provides Niv/14058, where we might slightly overpredict 
the observed emission. 

B.2. SMC stars 

AV 435 - 03 V((f )) (Fig. iBJOl The only discrepancies 
found during our analysis correspond to an overprediction of 
Heii/U6406, 6527, 6683. Both the Hei/Hen and the Nm/Niv 



Mas05: 02 If 
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ionization equilibrium suggest a hotter temperature than quoted 
by Mas05, = 46 kK. This temperature seems to be some- 
what cool for its spectral type 03 V assigned by Mas05 be- 
cause of NiV/14058 ^ N ni/14640, but quite consistent with our 
predictions for the derived wind-strength and nitrogen content, 
[N] =7.58 (cf. Figs. |7]and[S. 

AV 177 - 04 V((f)) (Fig. lBTTI The H/He analysis of this star 
produced similar parameters as found by Mas05. Owing to a 
high rotation, v sin i = 220 kms -1 , nitrogen lines are barely vis- 
ible in the spectrum. Weak traces of emission at N m/l/14634 - 
4640 - 4642 and N mU4510 - 4514 - 4518 together with weak 
absorption at NiV/16380 are fitted consistently at [N] = 7.78. 



NGC 346-355 - ON2 111(f) (Fig-ETJ). This st ar was consid- 



ered a s a standard for the 02 III(f*) category by Walborn et al 
( |2002|) , and later on updated to ON2 III(f ) by Iwalborn et a" 



(2004). As for Nl 1-031 (same type!) analyzed in paper II, we 
found problems to fit all N m/N iv/N v lines in parallel, but to a 
lesser extent. The basic difference is related to He 1/14471, which 
is not as clearly visible as for Nl 1-031. During our analysis, we 
considered two possible parameter sets: a cooler solution with 
T e s = 51 kK (red) and a hotter one with = 55 kK (black), us- 
ing either the N m/N iv or the N iv/N v ionization equilibrium. By 
inspection of He 1/14471, we note that both temperatures might 
be consistent with the very weak observed feature. For a sim- 
ilar nitrogen abundance, [N] = 7.98, we were able to fit either 
Niii/U4634 - 4640 - 4642, Nni/i/14510 - 4514 - 4518 and 
N iv/16380 for the cooler solution, or N iV/16380 and N V/L14603- 
4619 for the hotter one. 

Mas09, restricted to He 1/1447 1 as a primary temperature in- 
dicator, derived T e ff = 49.5 kK and logg = 3.9, which would 
agree with our cool solution, but is insufficient for the Niv/N v 
lines. The hotter solutio n is in better agre e ment with results 
from lBouret etal.1 (120031) and IWalborn et al.l (120041 ). who found 
r eff = 52.5 kK and logg = 4.0 fitting the Niv/Nv lines by 
mean s of cmfgen. In par ticular, we achieved a similar fit qual- 
ity as Bou ret et al.l d2003l) , for a similar nitrogen content. Bouret 
et al. stated that at r e jf ~ 55 kK (identical with our hotter solu- 
tion) their fit for He H/14686 would improve. Such an increase in 
temperature would also improve their fit of N v, which we are 
abl e to fit accurately. T he sa me stellar parameters as determined 
bv lBouret et al.l (120031) and IWalborn et al.l (2004) were derived 
by Heap et al. (2006j)using tlusty, mostly based on lines from 
highly ionized species, in particular N v and N iv. Unfortunately, 
they did not comment on Hei and Nm, but reassuringly they 
derived a nitrogen abundance very similar to ours, [N] = 7.92. 
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Fig. B.l. R136-040 - 02-3.5 V. Observed (green) and best fitting optical H/He and N spectrum. For the R136 O-stars, observed 
spectra for H a , H r , He 1/1/14387, 4471, 6678, Hen/14541, 6406, 6527, 6683, Nmvl/14510 - 4514 - 4518, and NiV/16380 taken from 
the STIS/CCD dataset. Remaining, lower quality spectra collected by FOS. Hg was not observed for this star. 
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B.6. LH 90:ST 2-22 - 03.5 III(f+). 
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g. B.12. NGC 346-355 - ON2 III(f ). Black: hotter solution (T eS = 55 kK), supported by N rv/N v lines. Blue: cooler solution (r eff 
51 kK), mostly supported by Nm (together with N rvA6380). Red: 'average' solution (r e (j = 53 kK) used in Sect. [7] 



